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We report a new value of electron g — 2, or a e , from 891 Feynman diagrams of order a 4 . 
The FORTRAN codes of 373 diagrams containing closed electron loops have been verified by 
at least two independent formulations. For the remaining 518 diagrams, which have no closed 



> 

■^j- \ lepton loop, verification by a second formulation is not yet attempted because of the enor- 

£-»« ' mous amount of additional work required. However, these integrals have structures that allow 

o" 

extensive cross-checking as well as detailed comparison with lower-order diagrams through the 

o 

renormalization procedure. No algebraic error has been uncovered for them. The numerical 



evaluation of the entire a 4 term by the integration routine VEGAS gives — 1.7283 (35) (a/ir) , 
where the uncertainty is obtained by careful examination of error estimates by VEGAS. This 
leads to a e = 1 159 652 175.86 (0.10) (0.26) (8.48) x 10~ 12 , where the uncertainties come from 
the a 4 term, the estimated uncertainty of a 5 term, and the inverse fine structure constant, 
a -1 = 137.036 000 3 (10), measured by atom interferometry combined with a frequency comb 
technique, respectively. The inverse fine structure constant a _1 (a e ) derived from the theory and 
the Seattle measurement of a e is 137.035 998 83 (51). 
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I. INTRODUCTION AND SUMMARY 



It was more than 18 years ago that an amazingly precise measurement of g — 2 values of 
the electron and positron in a Penning trap was reported Q]: 

aj(exp) = 1 159 652 188.4 (4.3) x 1(T 12 [3.7 ppb], 

a+(exp) = 1 159 652 187.9 (4.3) x 1(T 12 [3.7 ppb]. (1) 

Numerals 4.3 in parentheses represent uncertainties in the last digits of measurements, and 
3.7 ppb is the fractional precision (1 ppb = 1 x 1CU 9 ). At present a new measurement of 
a„ which is based on a cy.indrtca. cavity whose property has been studied aaa.ytically Q, 
is in progress. It is expected to improve the precision substantially [3J. Unlike the muon 
g — 2, which is sensitive to short distance physics, the electron g — 2 is rather insensitive 
to short distance physics. It is therefore calculable to very high precision from the QED 
or, more generally, the Standard Model. High precision measurements of a e thus enables us 
to subject the QED to a most stringent test. If disagreement is found between theory and 
experiment of a e , it would have a very profound impact on the validity of QED, or modern 
quantum field theory in general. This is why high precision measurement of a e is of such an 
importance. 

The QED contribution to a e can be written in general as 
a e (QED) = A 1 + A 2 {m e /m^) + A 2 (m e /m T ) + A 3 (m e /m^,m e /m T ). (2) 
Ai, A 2 , A 3 can be expanded as power series in a/n: 



^.^(2) + 4"(2) J + 4 B (2)' + .... -1.2.3. (3) 

Af\ A^\ and A^ are known analytically 0,0,0] whose numerical values are 
A? = 0.5, 

Aj 4) = -0.328 478 965 579 . . . , 

Af ] = 1.181 241 456 .... (4) 

Actually, the analytic result A^ is the culmination of a long sequence of analytic work 
that started around 1969 [l]. It took more than 25 years of hard work before it was finally 
completed by 0]. 



From the vantage of 1960 — 1970, analytic solution of the sixth-order muon g-2 diagrams 
containing an (external) light-by-light scattering subdiagram seemed to be extremely dif- 
ficult. Thus, one of the authors (T.K.) and Aldins developed a numerical approach based 
on the technique of Feynman-parametric integral [8|, calculating the integrand by hand. 
The subtle problem of renormalizing the light-by-light scattering subdiagram was bypassed 
taking advantage of gauge invariance. Meanwhile, Brodsky and Dufner were working on 
the same problem. We decided to join force and write the report together [9|. The result- 
ing 7-dimensional integral was evaluated by an adaptive-iterative Monte-Carlo routine jlol |. 
Sixth-order diagrams containing vacuum-polarization loops were also evaluated in a similar 
manner 

This method was then extended to general diagrams that require explicit renormalization 



12, 



131 ] . It consists of three steps: 

(1) Conversion of momentum space integrals into integrals over Feynman parameters. 
This step is fully analytic and carried out by an algebraic manipulation software, such as 
SCHOONSCHIP Q and FORM 0. 

(2) Renormalization by a subtraction term that subtracts divergence point- by-point 
throughout the domain of integration. The renormalization term is constructed in such 
a way that it is analytically factorisable into products of a divergent part of renormalization 
constant and a g-2 term of lower order. 

(3) Numerical integration by VEGAS [lfj], an improved version of adaptive-iterative 
Monte-Carlo integration routine jlo| . 

We applied this scheme to two different formulations, both in Feynman gauge. One is to 
apply it to each vertex individually and add them up. Another starts by combining a set of 
vertices into one with the help of the Ward-Takahashi identity 

^(p,g) = -£(p+|) + £(p-|), (5) 

where A M (p, q) is the sum of vertices obtained by inserting the external magnetic field in 
fermion lines of a self-energy diagram £(p). p±q/2 is outgoing (incoming) electron momen- 
tum. Differentiating both sides of (jSJ) with respect to q v one obtains 



A u (p,q) ~ -q^ 



dA^(p,q) 
dq v 



J g =o 



dp u 
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The magnetic projection of right-hand-side (RHS) provides the starting point of an inde- 
pendent formulation. This will be called Version A [17] . The approach starting from the 
left-hand-side will be called Version B. Although Version A requires additional algebraic 
work it results in fewer and more compact integrals than Version B and ensures significant 
economy of computing time. 

The entire a 3 term was coded in both Version A and Version B. By 1974 this led to a 
PI 

crude estimate of A\ 13]. To achieve a higher precision in this approach, however, very 
extensive computer work was required, which was not easily available at that time. It was 
only in 1995 that we were able to obtain a sufficiently precise result [18] 

Af ] (num.) = 1.181 259 (40). (7) 

Shortly afterwards this value was confirmed by the analytic result (f| to more than 5 digits. 

Although analytic technique developed for integrating A\ , in particular that of integra- 
tion by part, have been useful for analytic study of eighth-order term A\ ' further devel- 
opment seems to be needed for a complete analytic integration of A\ [19]. Only a small 
number of eighth-order diagrams without closed lepton loop have been integrated analyti- 
cally thus far 2f| . 

On the other hand, no theoretical difficulty has been encountered in_the Feynman- 
parametric approach to A{\ The renormalization scheme developed in 



12 1 for the sixth- 



order term could be readily extended to eighth-order. However, the enormous size of integrals 
of Group IV (see Sec. IVI|) and Group V (see Sec. IVII|) of eighth-order term and a large 
number of renormalization terms required presented a tremendous challenge to both alge- 
braic construction of integrands and their evaluation by numerical means. For these reasons 
numerical evaluation of these diagrams were carried out initially in Version A only, which 
required smaller codes than Version B. 

The first result of A® obtained by VEGAS 

4 8) = -1.434(138) (8) 

22j. This was a very preliminary result based on about 10 7 sam- 



was reported by 1990 21 



pling points per iteration and about 30 iterations. The estimated uncertainty of numerical 
evaluation was of the same order of magnitude as that of the experimental value 

In order to improve Ai beyond (jHJ) by an order of magnitude or more and move from 
a qualitative to quantitative calculation, it is necessary to evaluate the integrals using more 



Af 
Af 
Af 
Af 
Af 



than two orders of magnitude of sampling points. It was estimated that this would require 
computation of up to ten years on high speed computers then available. 

Because of this lengthy time scale and numerous requests for information on the status 
of calculation, preliminary values have been reported at several Conferences and Meetings 
and some of them have been printed in Proceedings and books: 

-1.557 (70) [23], 
-1.4092 (384), [24], 

-1.5098 (384), [25,26], ( 9 ) 

-1.7366 (60), [27], 
-1.7260 (50), [28,29]- 
Note that these values are not entirely independent. This is due to the fact that complete 
evaluation of each integral took typically three to six months of intense runs on fast computer. 
At any particular moment only a few of these integrals are being upgraded. The values in 
© are snapshots of this continuously evolving numerical work. 

(s) 

The difficulty encountered in improving numerical precision of A\ turned out to be far 
greater than that of Af\ It was partly due to the enormous size of integrands, which made 
accumulation of good sampling of integrand difficult. However, the most serious problem 
we have encountered is that our method is very sensitive to the round-off error (called digit- 
deficiency or d-d error) which is present in any computer calculation based on finite digit 
arithmetic, but affects our calculation particularly severely. (See Appendix B of Ref. 30j 
for details.) Because of these problems, some values given in (0) differ from each other 



24j differ by more than 2.6 s. d. Actually, 



considerably. For instance, the results of |25|] and 

error estimate of j^j generated by VEGAS was considerably smaller than that of |21j. 
However, unstable and erratic behavior of some integrals, in particular M40 discussed in 
detail in Sec. I VI II indicated that the error estimate generated by VEGAS could not be 

Qn 
was enlarged to that of |24|. 

More recently, we discovered an algebraic error in the integrands formulated in Version 
A Q while reevaluating 18 integrals of Group IV(d) by an independent formulation in 
Version B. Correction of this error is the cause of most of the jump from j^,!^ to j^f. 

The discovery of this error forced us to re-examine other diagrams, too. For this purpose, 
we have re-evaluated all diagrams of Groups IV(b) and IV(c) by an independent formulation 
in Version B. The result is that FORTRAN codes of all eighth-order g—2 diagrams containing 



closed electron loops (373 diagrams) have now been verified by more than one independent 
formulations (itIIsoI]. Remaining 518 diagrams belonging to Group V, which have no closed 
lepton loops, have thus far been formulated only in Version A, because the amount of work 
required for Version B would be several times larger than that of Version A. 

On the other hand, Group V diagrams are subject to very tight algebraic constraints. 
Namely, these diagrams have numerous renormalization terms which are constructed system- 
atically by a power-counting rule that enables us to carry out an extensive cross-checking 



22j . Having 



among themselves and with lower-order diagrams which are known exactly 
examined this procedure once again, we are convinced that they are completely free from 
algebraic error, although we found several typographical errors in the intermediate renormal- 



22 1 . These errors are corrected 



ization scheme described in Appendix C and Appendix D of 
in Appendix iBl of this paper. 

The r eally time-consuming task is to make sure that the result of numerical integration by 
VEGAS |1(t| can be trusted. This required a very good control of d-d error. (See Appendix 
B of Ref. |3j| for details.) This work has been carried out on a number of supercomputers 
over the period of more than 15 years. The latest results of this on-going effort are given 
in Eqs. (JMJ), (p?7j). (jSJ), (J52|) . and (f58j). Summing them up we obtain the total eighth-order 
term 



At 



(8) 



-1.7283(35), 



(10) 



where the uncertainty is determined by careful analysis of errors estimated by VEGAS. This 



25 



26). 



is an order of magnitude improvement over the previous value 

At the level of precision we are interested in we cannot totally ignore the tenth-order term 
A^°\ which is not known at present. In view of the gently increasing trend of coefficients of 
(a/rc) n as n increases from 1 to 4, however, it is unlikely that A] 10 ^ becomes suddenly very 
large. Following the argument of j3l|, let us assume, until proved otherwise, that 



4 10 ) = 0.0 (3.8) 



(11) 



In order to find out whether this is a justifiable assumption we have embarked on a very 



extensive project to evaluate A\ 



(10) 



32 
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The electron anomaly receives also small mass-dependent contributions 

Af(m e /m^ = 5.197 386 70 (27) x KT 7 , 
4 4) (m e /m T ) = 1.837 63 (60) x 10~ 9 , 
4 6) (m e /m M ) = -7.373 941 58 (28) x 10~ 6 , 
Af(m e /m T ) = -6.5819 (19) x 10" 8 , 
Af\m e /m^m e /m T ) = 0.190 95 (63) x KT 12 , (12) 

where the uncertainties are due to measurement uncertainties of m e /m ll and m e /m T only. 
Evaluation of A\ 8 ^(m e /m jU ), A^\m e /m T ), and Af\m e /m lJl ,m e /m T ) are not difficult since 
we already have their codes. But, numerically, they are totally insignificant at present. 

To evaluate the theoretical value of a e we must also include non-QED contributions of 
the Standard Model 

a e (hadron) = 1.671 (19) x 10" 12 , 
a e (weak) = 0.030 (1) x 10" 12 . (13) 

Although they are adequate for the moment, updates of these values would be desirable. 

In order to compare theory with experiment we need a value of a, and the most precise one 
'or that. Only recently such an a has become available by the progress of atom interferometry 
37| . Combined with the measurement of the cesium Dl line [38|, this leads to 

a-\h/M Cs ) = 137.036 000 3 (10) [7.4 ppb]. (14) 

Using this a and a e calculated from (JTUJ), (H2J), and (fTTj) we obtain 

a e (h/M Cs ) = 1 159 652 175.86 (0.10) (0.26) (8.48) x 10~ 12 , (15) 

where the first and second uncertainties are from (fTU|) and (fTTj) . and the third is from ()14j) . 
This leads to 

a e (exp) - a e (h/M Cs ) = 12.4 (4.3) (8.5) x 10~ 12 , (16) 

where first and second errors are from measurements of a e and a(h/Mc s ), respectively. 

It is noted in 0| that the uncertainty 7.4 ppb in ((T3J) may be brought down to 3.1 ppb 
when the nature of systematic error in the h/Mc s measurement is understood better. This 
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is very exciting: The discrepancy ()16|) between theory and experiment would then become 
2.2 s. d., which may provide the first serious challenge to the theory and measurement of a e . 
Of course further work is needed to find out whether this is a real discrepancy or something 
else. The result of new measurement of a e [| is anxiously waited for. 

An alternate and better test of QED is to compare a from (|14|) with a obtained from 
theory and measurement of a e : 

a-\a e ) = 137.035 998 834 (12) (31) (502) 

= 137.035 998 834 (503) [3.7 ppb], (17) 

where, in the first line, 12 and 31 are uncertainties in theoretical values (fTU|) and (jllj) and 
502 is that of experimental value This shows clearly that, if the measurement of a e is 
improved by an order of magnitude, the precision of a(a e ) will be enhanced to about 0.3 
ppb, putting it an order of magnitude ahead of a(h/M Cs ) of (JHJ). The current theoretical 
uncertainty in (a/ir^A^ is 0.10 x 10~ 12 , which is smaller than the uncertainty assumed for 
the (a /it) 5 term. This means that more accurate evaluation of will not improve the QED 
prediction significantly until a better estimate of A± is obtained. Complete evaluation of 
A^ 10 '', which has contributions from 12672 Feynman diagrams, each of which may occupy 
FORTRAN code up to two orders of magnitudes larger than those of eighth-order terms, 
clearly requires an enormous amount of work 

HQ- 
Appendix El describes how algebraic errors of Group V integrals are eliminated. Two-step 

renormalization scheme of Group V diagrams in Version A is described in Appendix [Bl 



II. CLASSIFICATION OF DIAGRAMS CONTRIBUTING TO Af } 

(8) 

A\ has contributions from 891 eighth-order vertex diagrams. Feynman integrals for 
these diagrams consist of twelve propagators integrated over 4 four- dimensional loop mo- 
menta. They may have subdiagrams of vacuum-polarization (v-p) type and/or light-by-light 
scattering (/-/) type. They fall into five gauge-invariant groups according to the type of 
proper vertex diagrams from which they are derived and the kind of closed electron loops 
they contain. The v-p diagrams found in A^ £1X6 clS follows: 
n 2 , which consists of one closed lepton loop of second-order. 
II4, which consists of three proper closed lepton loops of fourth-order. 
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n 4 (2), which consists of three lepton loops of type II4 whose internal photon line has a II2 
insertion. 

n 6 , which consists of 15 proper closed lepton loops of sixth-order. 

The l-l diagrams we need are: 
A 4 , which consists of six proper closed lepton loops of fourth-order, with four photon lines 
attached to them. 

(2 s ) 

, which consists of 60 diagrams in which lepton lines and vertices of A 4 are modified by 
second-order radiative corrections. 

We are now ready to classify the diagrams into five (gauge- invariant) groups: 
Group I. Second-order electron vertex diagrams containing lepton v-p loops n 2 , n 4 , IT 4 ( 2 ) 
and/or Il 6 . This group consists of 25 diagrams. 

Group II. Fourth-order proper vertex diagrams containing lepton v-p loops H2 and/or IT 4 . 
This group consists of 54 diagrams. 

Group III. Sixth-order proper vertex diagrams containing a v-p loop 1I2. This group 
consists of 150 diagrams. 

(2) 

Group IV. Vertex diagrams containing a light-by-light scattering subdiagram A y i J or A 4 
with further radiative corrections of various kinds. This group has 144 diagrams, which can 
be divided further into subgroups IV(a), IV(b), IV(c), and IV(d). 
Group V. Vertex diagrams with no closed electron loop. It consists of 518 diagrams. 

Diagrams of Groups I, II, III and IV(a) are relatively easy to deal with since their struc- 
tures are closely related to the second-, fourth-, and sixth-order diagrams, which are known 
exactly. They have been checked by several independent calculations j^J. 

In order to deal with the remaining diagrams in general, we need a full power of algebraic 



formulation originally developed for sixth-order diagrams 



order diagrams in 



12j and extended to the eighth- 



2l\ . But, for diagrams of Groups IV(b), IV(c), IV(d), as well as Group 



V, the Version B codes are so huge that we chose initially to work only with Version A. 
Since we fixed the error in Group IV(d) with the help of Version B calculation, however, we 
decided to evaluate Groups IV(b) and IV(c) in Version B, too. The numerical results are in 
excellent agreement with Version A jscj. 

This leaves Group V as the only group which has been evaluated in Version A only. Until 
it is evaluated in Version B or other methods, we have to rely on the very tight algebraic 
construction as the guarantor of the FORTRAN code of Group V. 
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III. GROUP I DIAGRAMS 



Group I diagrams can be classified further into four gauge-invariant subgroups: 

Subgroup 1(a). Diagrams obtained by inserting three n 2 's in a second-order vertex. This 
subgroup consists of one Feynman diagram. See Fig. Ufa). 

Subgroup 1(b). Diagrams obtained by inserting a II2 and a II4 in a second-order vertex. 
Six Feynman diagrams belong to this subgroup. See Fig. Gib). 

Subgroup 1(c). Diagrams containing n 4 ( 2 ). There are three Feynman diagrams that belong 
to this subgroup. See Fig. |21 

Subgroup 1(d). Diagrams obtained by insertion of Il 6 in a second-order electron vertex. 
Fifteen Feynman diagrams belong to this subgroup. Eight are shown in Fig. |21 Diagrams 
a, c, d, e, / and the time-reversed diagram of e have charge-conjugated counterparts. 

The evaluation of subgroups 1(a) and 1(b) is greatly facilitated by the analytic formulas 

available for the second- and fourth-order spectral representations of the renormalized photon 

I I 

propagators [29j . The contribution to a e from the diagram obtained by sequential insertion 
of m k-th order electron and n Z-th order electron v-p loops into a second-order electron 
vertex is reduced to a simple formula 



« = / dy(l- y) 
'0 



ds 



Pk{s) 



y m e 



m,. 



1 

dt — 
1 



t 2 y 2 



(18) 



where pk is the fc-th order photon spectral function. Exact p 2 and p 4 can be found in Ref. 



21 
^0 



39 
ill- 



. An exact spectral function for 114(2) and an approximate one for Hq are also available 



The contribution of diagrams of Fig. Qcan be obtained by choosing (k = 2, m = 3, n = 
0), (k = 2, m — 2, 1 = 2, n — 1), (k = 2, m — 1, 1 = 2, n = 2). The latest numerical values 
obtained by evaluating these integrals using VEGAS are listed in Table HI where the 
number of sampling points per iteration and the number of iterations are also listed. 

Note that these diagrams need no additional renormalization. Thus the renormalized 
amplitude 02,^2:3 is given by 



°2,p2:3 — M 2 , p 2:3- 



(19) 
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i(b) 



FIG. 1: (a) Diagrams contributing to subgroup 1(a). (b) Diagrams contributing to subgroup 1(b). 
Solid horizontal lines represent the electron in external magnetic field. Numerals "2", "4" within 
solid circles refer to the proper renormalized v-p diagrams II2 and II4 , respectively. One and six 
Feynman diagrams contribute to 1(a) and 1(b), respectively. 



From Table HI we obtain the contribution of subgroup 1(a) 



0.000 876 911 (40), 



(20) 



which is 5 times more precise than the earlier result 
the analytic result 



42j , and is in excellent agreement with 



a f) A anal. 



0.000 876 865 



(21) 



The contributions of Fig. ^b) for (li, I2) = (e, e) can be written down in a similar fashion. 
The most recent result of numerical integration by VEGAS is listed in the second row of 
Table □ 



i {8) 
h(b) 



0.015 325 6 (6). 



(22) 



M 

FIG. 2: Diagrams contributing to subgroup 1(c). (h^h) = (e, e). See FIG. 1 for notation. 

In evaluating the contribution to a e from the 9 Feynman diagrams of subgroup 1(c) shown 
in FIG. [21 our approach is to make use of the parametric integral representation of the v-p 
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TABLE I: Contributions of diagrams of Figs, ^a), ^b), and [3 np is the number of Feynman 
diagrams represented by the integral. Evaluation was carried out on a workstations in 1997. 



Integral 


iif 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


-^2,P2:3 


1 


0.876 911 (40) x 10~ 3 


1 x 10 7 


55 


M2,P2,P4 


6 


0.015 325 6 (6) 


1 x 10 7 


55 


AM 2i p 4a (P2) 


1 


0.011 403 1 (8) 


1 x 10 7 


64 


AM 2j Pib(P2) 


2 


0.001 716 8 (4) 


1 x 10 7 


55 


AM 2jP6a 


2 


0.044 448 (3) 


1 x 10 8 


70 


AM 2iP6b 


1 


0.028 594 (3) 


1 x 10 8 


70 


AM 2 , P6c 


2 


-0.038 374 (2) 


1 x 10 8 


70 


AM 2jP6d 


2 


-0.027 507 (2) 


1 x 10 8 


70 


AM 2jP6e 


4 


0.179 335 (3) 


4 x 10 8 


70 


AM 2)P6f 


2 


-0.062 004 (2) 


1 x 10 8 


70 


AM 2jP6s 


1 


0.038 879 (2) 


1 x 10 8 


70 


AM 2jP6/l 


1 


0.023 675 (2) 


1 x 10 8 


70 



term LT4( 2 ). Following the two-step renormalization procedure, these contributions can be 
written in the form |44( 

G2,P4(P2) = AM 2j p 4a (p 2 ) + AM 2) p 4fc (p 2 ) — 2A£? 2) p 2 M 2i p 2 , (23) 

where each term is finite integral obtained by the K5 renormalization procedure described in 
Ref. 01 • The suffix P2 stands for the second-order v-p diagram LT 2 , P4 for the fourth- 
order v-p diagram LT4, while P4(P2) represents the diagram n 4 ( 2 ). P4 receives contributions 
from P^a (vertex correction) and P46 (two lepton self-energy insertions), P4 = P± a + P4&. 

The results of numerical evaluation of (|2*3*|). obtained by VEGAS, are listed in Table |TJ 
Note that the value listed in column 3 for AM 2)P m P2 ) includes the factor n F . Numerical 

12 



TABLE II: Auxiliary integrals for Group I. Some integrals are known exactly. Remaining integrals 
are obtained numerically by VEGAS. 

Integral Value (Error) Integral Value (Error) 

M 2 ,p 2 0.015 687 421 •• • M 2jP2 . -0.012 702 383 ■■ ■ 

AAf 2 ,p 4 0.076 401 785 ■• • 

AB 2 0.75 AB 2 , P2 0.063 399 266 ■• • 

AL 4 0.465 024 (17) AB 4 -0.437 094 (21) 

A5m 4 1.906 340 (21) 

values of lower-order Feynman integrals, in terms of which the residual renormalization 
terms are expressed, are given in Table |HI From these Tables we obtain 

af} = 0.011 130 8 (9). (24) 

We obtained an independent check of (J24j) using a two-dimensional integral form of a 3 
spectral function for LT 4 ( 2 ) of Fig. El which was derived from the QCD spectral function 
obtained in [40] . Numerical integration using this spectral function gives 



i (8) 

l Ic 



0.011 131 2 (12), (25) 



for 10 million sampling points iterated 20 times in quadruple precision. This is in agreement 
with (J24p to the fifth decimal points althou gh t heir approaches are quite different. 
The new result ()24|) confirms the old result |22| but with a much higher precision. 

The contribution to a e from 15 diagrams of subgroup 1(d) (see Fig. EJ) can be written as 

a>2,P6i = AAf 2) P6i + residual renormalization terms , (i = a, . . . , h). (26) 



Divergence-free integrals AM 2j p 6i are defined by (4.13) of Ref. 2l|. Their numerical values 
(summed over the diagrams related by time-reversal and charge-conjugation symmetries) 
are evaluated numerically by VEGAS and listed in the third column of Table HI 

Summing up the contributions of diagrams P6a to P6h of Fig. 01 we obtain the following 
expression: 

h 

a Hd) = AM 2,p 6 * - 4A J B 2 AM 2i p 4 

i=a 

+ 5(AB 2 ) 2 M 2j p 2 - 2(AL 4 + A5 4 )M 2i p 2 

- 2A5m 4 M 2 ,p 2 *, (27) 
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P6a 


P6b 


P6c 


P6d 


P6e 


P6f 


P6g 


P6h 



FIG. 3: Eighth-order vertices of subgroup 1(d) obtained by insertion of sixth-order (single electron 
loop) v-p diagram Ii6 in a second-order electron vertex. 



where the value of AM 2> pq{ m column 3 of Table |l] is the sum of rip vertex diagrams and 

AB 2 = A'B 2 + A'L 2 = -, 

AM 2iP4 = AM 2>P4a + 2AM 2jP4fe , 

AL 4 = AL 4x + 2AL Ac + AL Al + 2AL 4s , 

AB 4 = AB Aa + AB 4b , 

A8m A = A5m Aa + A5m A b- (28) 

The quantities listed in (|2*H|) are defined in Ref. [21]. Their numerical values are listed in 
Table |H] The 1998 results of numerical integration of AM 2) p6i are listed in Table HI From 
the numerical values in Tables HI and HP we obtain 

af^ = 0.049 516 1 (63). (29) 

We obtained an independent check of (|29jl using the Pade approximant of the sixth-order 
photon spectral function 

af {d) (Pade) = 0.049 519 5 (38), (30) 

which is in excellent agreement with (|29|). Calculation was carried out in real*16 for 10 
million sampling points iterated 100 times. 

Collecting the results (|2*T|) . (|2*2*j) . (|2*4^l and (|2*9*)l . we find the best value of the contribution 
to the electron anomaly from the 25 diagrams of group I to be 

Af ] [I] = 0.076 849 (7). (31) 
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IV. GROUP II DIAGRAMS 



Diagrams of this group are generated by inserting II2 and II4 in the photon lines of fourth- 
order vertex diagrams. Use of analytic expressions for the second- and fourth-order spectral 
functions for the photon propagators and time-reversal symmetry cuts down the number of 
independent integrals in Version A from 54 to 8. 

The contribution to a e arising from the set of vertex diagrams represented by the "self- 
energy" diagrams of Fig. 0] can be written in the form 

a 4,p a — AM 4j p Q + residual renormalization terms, (32) 

where AM^p a are finite integrals obtained in the intermediate step of two-step renormaliza- 
tion 22]. Their numerical values, obtained by VEGAS are listed in Table ITTT1 The values of 
auxiliary integrals needed to calculate the total contribution of group II diagrams are given 
in Tables El and HY1 




FIG. 4: Eighth-order diagrams obtained from the fourth-order vertex diagrams by inserting 
vacuum-polarization loops II2 and II4. 

Summing the contributions of diagrams of the first, second, and third rows of Fig. 01 one 
obtains (see Ref. (2^] for notation) 

a 4iP4 = AM 4a)P4 + AM 4bPl ,. A + AM 46i p :4 

- AB 2 M 2jP4 - AB 2iP4 M 2 , (33) 



15 



TABLE III: Contributions of diagrams of Fig. 0] np is the number of Feynman diagrams rep- 
resented by the integral. Suffixes PO and PI' are remnants of notations in 22] and refer to two 



virtual photons in which v-p insertions are made. Numerical evaluation was carried out on a 
workstations in 2001. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


AM 4a ,P4 


18 


0.131 298 (8) 


2 x 10 9 


100 


AM 4f)i p 0: 4 










+AM 4b> Pv. A 


18 


-0.420 295 (8) 


4 x 10 8 


100 




3 


-0.004 351 (2) 


1 x 10 8 


70 


AMi&,Pl':2,PQ:2 


3 


-0.022 332 (1) 


1 x 10 8 


70 




6 


0.007 875 (8) 


2 x 10 8 


70 


AM4ft j po:2:2 










+ AM 4 b j p 1 / : 2:2 


6 


-0.065 458 (4) 


1 x 10 8 


70 



0-4,P2,P2 — AM4a,P2,P2 + AM 46 pi' : 2PQ : 2 

- A£ 2 ,P2 M 2 ,p 2 , (34) 

Q-4,P2:2 = AM 4a] P2:2 + ^M4i, i Pl l : 2:2 + ^^4d,P0:2:2 

- AB 2 M 2tP2 .. 2 - Afi 2 ,P2:2 M 2 , (35) 

respectively, where M 2j p 4 is equal to AM 2j p 4 — 2Afi 2 A^2,P2- M^.. corresponds to the left- 
most diagrams of Fig. 0] and M^,.. ones on the right. AM 4ai ... is the sum of np vertex 
diagrams. On the other hand, the auxiliary integrals listed in Tables |H] and II VI do not 
include multiplicity. 

Substituting the data from Tables ILT1 and IIHI into (jSSJ)? (jMj) , and (|33|) we obtain 

a 4P4 = - 0.420 483 (12), 
a4,P2,p 2 = - 0.027 677 ( 2), 
a4,P2:2 = - 0.073 452 ( 9). (36) 

Note that values listed in column 3 of Table ILTT1 include the factor n F . 
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TABLE IV: Auxiliary integrals for Group II. Some integrals are known exactly. Remaining inte- 
grals are obtained by VEGAS integration, with total sampling points of order 10 11 . 



Integral 


Value (Error) 


Integral 


Value (Error) 


M 2 


0.5 


M 2 ,P4 


0.052 870 652 •• • 


M2,P2:2 


0.002 558 524 •• • 






Ai? 2 ,P4 


0.183 666 5 (18) 




0.027 902 3 (4) 



The sum of contributions of 54 diagrams of group II is 
a { % = - 0.521 613 (15). (37) 

V. GROUP III DIAGRAMS 

Diagrams belonging to this group are generated by inserting a second-order vacuum- 
polarization loop n 2 in the photon lines of sixth-order electron vertex diagrams of the three- 
photon-exchange type. Time-reversal invariance and use of the function p 2 (see ()18j) ) for the 
photon spectral function reduce the number of independent integrals in Version A from 150 
to 8. Some of these integrals are represented by the "self-energy" diagrams of Fig. |S| 

Let Mq 0) p2 be the magnetic moment projection in Version A of the set of 150 diagrams 
generated from a self-energy diagram a (=A through H) of Fig. El by insertion of LT2 and an 
external vertex. The renormalized contribution due to the group III diagrams can then be 
written as 

H 

'■III 



afrr — Ct6a,P2, (38) 



where 

a6«,P2 = AM 6q , i p2 + residual renormalization terms. (39) 

AM 6q , 5 p2 is UV- and IR-finite: all divergences of M 6Qj p 2 have been projected out by and 
1r operations. (See Ref. 13 •) 

The latest numerical values of Group III integrals are summarized in Table IV! Numerical 
values of auxiliary integrals needed in the renormalization scheme are listed in Tables |Hj IIVI 
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6H.P2 



FIG. 5: Typical eighth-order diagrams obtained by insertion of a vacuum-polarization loop H2 in 
electron diagrams of the three-photon-exchange type. Altogether there are 150 diagrams of this 
type. 



and ED For comparison, the results of old calculation |42j carried out in double precision are 
listed in the last column of Table This is to examine the effect of digit- deficiency error. 
In this case the effect is relatively mild because the introduction of a vacuum-polarization 
loop tends to make the integrand less sensitive to the singularity. 

When summed over all the diagrams of group III, the UV- and IR-divergent pieces cancel 
out and the total contribution to a e can be written as a sum of finite pieces (see Ref. 
for notation): 



H 



i (8) 
l III 



^AM f 



6a,P2 



a=A 



- 3AB 2 , P2 AM 4 - 3AB 2 AM 4 ,p 2 

+ (M 2 , iP2 [J] - M 2 , iP2 )A<5m 4 + (M 2 *[I] - M 2 *)A8m i)P2 

- M 2 , P2 [AB 4 + 2AL 4 - 2(AB 2 ) 2 } 

- M 2 (AB 4tP2 + 2AL 4 , P2 - 4AB 2 AB 2tP2 ), 



(40) 



where AM 6a p 2 listed in Table IS is the sum of n P vertex diagrams. M 2 * is the magnetic 
moment obtained from M 2 by insertion of a two-point vertex in the electron lines. M 2 * [I] is 
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TABLE V: Contributions of diagrams of Fig. \5\ rap is the number of Feynman diagrams rep- 
resented by the integral. Integrals were evaluated in double precision except for a small part of 
AMqg,P2 which was evaluated in quadruple precision. 



integral 


rap 


Value (Error) 


Sampling per 


JNo. of 


Data from 






including rap 


iteration 


iterations 


Ref. [42] 


AM 6 A,P2 


5 


-0.440 307 (56) 


1 x 10 8 


97 


-0.440 3 (3) 


AAf 6BiP2 


5 


0.689 660 (84) 


2 x 10 8 


100 


0.690 7 (7) 


AM 6C ,P2 


5 


0.638 503 (80) 


2 x 10 8 


100 


0.638 2 (7) 


AM 6AP2 


10 


0.440 999 (77) 


2 x 10 8 


100 


0.442 (7) 


AM 6B , P2 


5 


0.431 541 (55) 


1 x 10 8 


100 


0.432 1 (4) 


AM 6F<P2 


5 


0.442 163 (83) 


2 x 10 8 


100 


0.441 9 (8) 


AM 6G) p 2 


10 


1.140 811 (76) 






1.140 4 (10) 


d-part 




1.054 635 (60) 


2 x 10 8 


100 




q-part 




0.086 176 (47) 


2 x 10 7 


100 




AM 6 p iP2 


5 


-0.797 941 (83) 


6 x 10 8 


100 


-0.797 5 (12) 



TABLE VI: Auxiliary integrals for Group III. Some integrals are known exactly. Remaining 
integrals are obtained by VEGAS integration, with total sampling points of order 10 11 . 



Integral Value (Error) Integral Value (Error) 

M 2 * 1.0 M 2 *[I] -1.0 

M 2 *,p 2 0.044 077 4 (3) M 2 *, P2 [I] 0.010 255 3 (11) 

AM 4 0.030 833 612 •• • AM 4)P2 -0.106 707 082 ■• • 

A<5m 4 1.906 340 (21) /\5m AyP2 0.679 769 (15) 



AB 4)P2 + 2AL 4iP2 0.085 865 (30) 



the part of M 2 * which may be IR-divergent according to the IR power counting. An analytic 
evaluation of M 2 * [I] turns out to be finite and is listed in Table IV11 
Plugging the values listed in Tables ILTl II VI and I VII in (|4"0|). we obtain 



a$j = 1.417 722 (214). (41) 
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VI. GROUP IV DIAGRAMS 



Diagrams of this group can be divided into four subgroups: IV(a), IV(b), IV(c), and 
IV(d). Each subgroup consists of two equivalent sets of diagrams related by charge conjuga- 
tion (reversal of the direction of momentum flow in the loop of the light-by- light scattering 
subdiagram). Diagrams of subgroups IV(a), IV(b), and IV(c) are obtained by modifying 
the sixth-order diagram Mqh which contains the light-by-light scattering subdiagram A 4 , 
one of whose external photon line represents the magnetic field. The magnetic moment 
contribution Mqh is known analytically 45^ whose numerical value is 

M 6LL = 0.371 005 292.... (42) 

Subgroup IV(a). Diagrams obtained by inserting a second-order vacuum-polarization 
loop H2 in Mqll. This subgroup consists of 18 vertex diagrams. In Version A they are 
represented by the self-energy-like diagrams shown in Fig. El Insertion of the magnetic 
vertex in these figures generates diagrams of Version B. Since M^n is known analytically 
and since insertion of H2 is straightforward, this subgroup is known exactly. 




FIG. 6: Electron self-energy- like diagrams representing the external- vertex-summed integrals of 
subgroup IV(a). (li,h) = (e,e). 

Subgroup IV(b). Diagrams containing sixth-order light-by- light scattering subdiagrams 
A 4 ( 2 ). Altogether, there are 60 diagrams of this type. Charge-conjugation and time-reversal 
symmetries and summation over external vertex insertions reduce the number of independent 
integrals to 4 in Version A, represented by the self-energy- like diagrams LLA, LLB, LLC 
and LLD of Fig. They can be generated from a common template by permutation of 
tensor indices of virtual photons attached to the A 4 loop. Similarly, all diagrams of Version 
B can be generated from a common template. 
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LLA LLB LLC LLD 




FIG. 7: Electron self-energy-like diagrams representing (external-vertex-summed) integrals of 
subgroup IV(b), IV (c), and IV(d). 

Subgroup IV(c). Diagrams obtained by including second-order radiative corrections on 
the open electron line of M^^. There are 48 diagrams that belong to this subgroup. Sum- 
mation over external vertex insertions and use of the interrelations available due to charge- 
conjugation and time-reversal symmetries leave five independent integrals in Version A, 
represented by the self-energy-like diagrams LLE, LLF, LLG, LLH and LLI of Fig. [7| They 
can be generated from a common template by permutation of tensor indices of virtual pho- 
tons attached to the open electron line. Similarly, all diagrams of Version B can be generated 
from a common template. 

Subgroup IV(d). Diagrams generated by inserting A4 internally in fourth-order vertex 
diagrams. Diagrams of this type appear for the first time in the eighth order. Charge- 
conjugation invariance and summation over the external vertex insertion with the help of the 
Ward-Takahashi identity lead us to three independent integrals in Version A, represented by 
the diagrams LLJ, LLK and LLL of Fig. [7| They can be generated from a common template 
by permutation of tensor indices of virtual photons attached to the A4 loop. Similarly, all 
diagrams of Version B can be generated from a common template. 

In subgroups IV(a), IV(b), and IV(c) UV-divergences arising from the light-by-light scat- 
tering subdiagram A 4 , or more explicitly W a ^(q, k h kj, fcj), 

can be taken care of by making 
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use of the identity: 



(43) 



which follows from the Ward-Takahashi identity. Namely, no explicit UV renormalization 
is needed if one uses the RHS of (}4*3*|) instead of LHS and the fact that S(p) of (jUJ) vanishes 
by Furry's theorem. On the other hand, S(p) is nonzero for subgroup IV(d) and the UV- 
divergence associated with the light-by-light scattering subdiagram A4 must be regularized, 
e. g., by dimensional regularization. For these diagrams it is necessary to carry out explicit 
renormalization of A 4 as well as that of the two sixth-order vertex subdiagrams containing A4. 
See fo r a detailed discussion of renormalization based on a combination of dimensional 
regularization and Pauli-Villars regularization. 

Many diagrams of Groups IV(b) and IV(c) have UV divergences associated with the 
second-order vertex or self-energy parts. Their renormalization is achieved by a subtraction 
term which can be factorized analytically into a product of UV-divergent part of the renor- 
malization constant and the sixth-order amplitude M§ll (in Version A) or its component (in 
Version B). This enables us to verify the analytic structure of many terms of the integrand. 
A similar argument applies to diagrams with IR divergence. Since different diagrams gen- 
erated by the same template have different divergences, this procedure explores the validity 
of template itself exhaustively. 

On the other hand, Group IV(d) diagrams LLJ, LLK, and LLL have no UV-divergent 
subdiagram of second-order. Thus they could not be verified by comparison with known 
lower-order structure. Indeed, this is why the errors in 18 integrals of Group IV(d) remained 
undetected until the integrand was constructed in Version B The error was in the 

template of Group IV(d) where some non-divergent terms were left out inadvertently. 

The argument presented above suggests that the validity of codes for Groups IV(b) and 
IV(c) can be established within Version A. However, as a further check, we have also evalu- 
ated them in Version B. In the following let us consider Version A and Version B separately 
since renormalization is handled slightly differently in two cases. 
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A. Version A 



simple. This is because Mqll 
, and insertion of the vacuum 



The calculation of group IV(a) contribution is particular 
has been fully tested by comparison with the analytic result 
polarization term is straightforward. As a consequence the integral M 6 ££ p 2 is finite, which 
means 



G6LL,P2 



M, 



6LL,P2 



AM 



6LL,P2- 



(44) 



The result of numerical calculation is listed in Table IVII 

Let us denote magnetic projections of subgroups IV(b) and IV(c) as M 8LLa where a 



A, ..,1. Relating the IR- and UV-divergent M 8LLa to the finite, numerica. 
AMsLLa defined by the procedure of two-step renormalization of Ref. 
the contributions of the diagrams of subgroups IV(b) and IV(c) as 

D 



and 



i {8) 



l {8) 
IV (c) 



J2&M 8 LLa - 3AB 2 M 6LL , 



ly calculable piece 
4611 . one can write 



(45) 



a=A 



^AM 8LiQ - 2AB 2 M 6LL , 



a=E 



(46) 

and AM 8LLa 



where the six-order light-by-light contribution to the anomaly is given by 
listed in Table lyTIl represents the sum of rip diagrams. 

Numerical integration of all terms contributing to a\y has been carried out using VEGAS 
jlfil ]. The latest results for Groups IV(b) and IV(c) are listed in Table IVTT1 The result for 
Group IV(d) had been handled separately In general, the major difficulty in dealing 
with the diagrams of Groups IV(b) and IV(c) arises from the enormous size of integrands 
(up to 5000 terms and 240 kilobytes of FORTRAN source code per integral) and the large 
number of integration variables (up to 10). 

Unlike A 2 (m^/rrie) for muon g — 2, which has a singular surface just outside of the 
integration domain (unit cube), the singularity of the mass-independent A\ is far away 
from the domain of integration. This makes evaluation of electron g — 2 less sensitive to the 
digit- deficiency problem. Nevertheless, our point-wise renormalization scheme is susceptible 
to the d-d problem. 
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TABLE VII: Contributions of dia gram s of Fig. Eland Fig. Q excluding LLJ, LLK, and LLL which 
were evaluated separately in Ref. 



171 ] . np is the number of Feynman diagrams represented by the 
integral. Some integrals are split into two parts: d-part is evaluated in real*8 and q-part is evaluated 
in real*16. a-part refers to the adjustable precision method developed by 4^|- The superscript * 
indicates that indicated contributions were obtained by extrapolation from calculations in which 



301 ] for details. 



the edges of integration domain were chopped off by l.d-10. See Appendix B of Ref. 
Numerical work was carried out on SP3, velocity cluster, SP2, Condor cluster, and a workstations 
over several years. The table lists only the latest of results obtained by various means. 



Integral 


np 


Value (Error) 


Sampling 


per 


No. of 






including np 


iteration 


iterations 


AMfiTT po 


6 


0.598 948 (137) 


4 x 


10 s 


100 


AMsr j A 


10 


0.171 527 (292) 








d-part 




0.148 729 (268) 


2 x 


10 8 


100 


n-rtfivf 

\J LJVAj 1 Li 




022 798 (116) 


2 x 


10 5 


380 


AM 8LLB 


20 


- 1.296 116 (481) 


2 x 


10 8 


99 


AM 8LLC 


20 


2.840 029 (513) 








d-part 




2.705 399 (345) 


1 X 


10 9 


60 


q-partl 




0.132 081 (375) 


4 x 


10 6 


63 


q-part2 




0.002 549 ( 53) 


1 X 


10 6 


100 


AM 8LLD 


10 


- 0.058 890 (325) 


2 x 


10 8 


100 


AM 8LLE 


6 


- 0.265 962 (201) 


1 X 


10 8 


80 


AM 8 LLF 


12 


- 0.764 679 (413) 








d-part 




- 0.618 362 (346) 


2 x 


10 8 


80 


q-partl 




-0.142 801 (221) 


2 x 


10 7 


120 


q-part2 




- 0.003 516 ( 36) 


1 X 


10 6 


400 


AMsllg 


12 


- 0.550 968 (670) 


4 x 


10 8 


80 


AM SLLH 


6 


0.190 770 (682) 








d-part 




0.069 275 (379) 


2 x 


10 8 


100 


q-part 




0.121 495 (566) 


2 x 


10 7 


100 


AM 8LLI 


12 


0.808 638 (482) 


4 x 


10 8 


80 
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In most cases the first step is to make the integrand smoother by stretching (see Appendix 
B of Ref. jscj), which is repeated several times until the integrand behaves more gently. 

Although chopping (see Appendix B of Ref. j^) was handy to obtain a rough estimate 
quickly, we had to abandon it in the end because extrapolation to 5 = turned out to be 
too unreliable in order to reach the desired precision. 

Most integrals were then evaluated by splitting them into two parts, one evaluated in 
real*8 and the other in real*16. In some cases, however, even the part evaluated in real*16 
suffered from severe d-d problem, preventing us from collecting large enough samplings for 
high statistics. Analyzing this problem closely, we found that it is possible to evaluate 
these integrals by the following procedure: First try several iterations with a positive rescale 
parameter (3 (typically (3 = 0.5) until VEGAS begins to show sign of blowing up due to the 
d-d problem. Then freeze [3 to (see Appendix B of Ref. |2oJ|). This may solve the problem 
in most cases. If not, try several iterations and see how rapidly the calculation runs into the 
d-d problem. It turns out that it takes place very early if we chose too many sampling points 
Ng per iteration. This is because choosing large N$ increases the chance of hitting random 
numbers too close to the singularity within one iteration. As a consequence the d-d problem 
is likely to dominate each iteration and makes it very difficult to collect large enough number 
of good samplings. A better strategy would be to reduce the size of Ns to a moderate value 
and, instead, increase the number of iterations Nj substantially. This is acceptable since, 
for [3 = which means that the distribution function p is no longer changed from iteration 
to iteration, the final error generated by VEGAS depends only on the product NsNi. 

This strategy has been applied in particular to the diagrams LLA and LLC, as is seen 
from Table IVTT1 Entries in Table IVTT1 are only the best of results obtained by various methods 
discussed above. They are consistent with each other despite their diverse approaches. 

One obtains from Table IVTT1 the Version A contributions of subgroups IV(b) and IV(c): 



a 



(8) 

'IV (b) 



0.821 79 ( 83) 



a 



(8) 

>IV(c) 



1.138 71 (117). 



(47) 



B. Version B 



In Version B the magnetic moment projection is evaluated for each vertex diagram on the 
LHS of ©. It is convenient to denote these diagrams in terms of self-energy-like diagrams 
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of Fig. [3 by attaching suffix % to indicate the lepton line in which an external magnetic field 
vertex is inserted. For instance, we obtain vertex diagrams LLA1, LLA2, LLA5 from the 
diagram LLA. 

We will not discuss subgroup IV (a) here since its Version A has already been fully tested. 
For subgroup IV(b) we find 

D 5 

a fv(b) = EE AM ^« " 4A£? 2M 6LL , (48) 

a=A i=l 

instead of ()45j) . Note that the last term of (|48|) is different from that of (|45p. This is not an 
error. It arises from difference in the definition of AM terms. Similarly, for subgroup IV(c) 
we obtain 

/ 3 

4? (C )= EE M ^« " 2AB ^eu- (49) 

a=E i=l 

Numerical evaluation of AM 8 n ai are carried out with 10 9 sampling points per iteration. 
The results are listed in Table IVIIII Precision of these calculations is somewhat higher 
than that of Version A. See the last column of Table IVTTTl for comparison of two Versions. 
Although the values of AM 8LL a and AM 8LLC in two versions differ by several s. d., this is 
due to the influence of d-d problem. The general agreement between Version A and Version 
B leaves no doubt that their codes are free from any error. 

One obtains from Table IVTTTl the values of afy^ and afy,, 

af v(b) = 0.822 57 (30), 

af v{c) = -1.138 93 (37), (50) 
which are consistent with those given in (|47|). 

C. Total contribution of Group IV 

The contribution of subgroup IV(a) is listed only in ()47|) since it was not evaluated in 
Version B. Previously only subgroup IV(d) was evaluated in both versions with comparable 
statistical weights [17 1. Now that subgroups IV(b) and IV(c) have been evaluated in two 
independent ways, we may combine their results statistically. This leads to our best estimate 
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TABLE VIII: Contribution of Group IV(b) and Group IV(c) diagrams of Fig. [7| evaluated in 
Version B. Double precision is used for all calculations which are carried out on RSCC at RIKEN. 
Finite renormalization terms AB%M%i,i,i^i = 1,2,3, are needed for LLA and LLC, respectively, in 
order to compare them with the calculations in Version A. M 6 ^^( 2 +3) = ^6LL2 + Mblz,3 is obtained 
subtracting Mqlli from the known value of Mqll(= Mqlli + A#6ll2 + Mqllz) given in l|l2*|). 



Integral 




Value (Error) 


Sampling 


per 


No. of 




Difference 






including np 


iteration 


iterations 




Ver. 


A - 


Ver.B 




10 


0.171 789 (131) 










-0. 


nnn 
UUU 




[oZV ) 


7] A AM 8 LLAi 




-0.096 030 (128) 


1 


X 


10 M 


300 










-AB 2 M 6LL2 




0.267 819 ( 27) 


1 


X 


10 9 


100 










AM 8UB 


20 


-1.296 748 (161) 


1 


X 


10 9 


350 


+0 


.000 


631 


(507) 


AM 8LLC 


20 


2.840 956 (196) 










-0 


.000 


927 


(549) 


Ei=l VC^MsLLCi 




3.387 028 (194) 


1 


X 


10 9 


516 










-AB 2 M 6ii ( 1+3 ) 




-0.546 072 ( 27) 


1 


X 


10 9 












AM 8LLD 


10 


-0.058 661 ( 72) 


1 


X 


10 9 


420 


-0 


.000 


229 


(332) 


AM 8LLS 


6 


-0.265 908 ( 84) 


1 


X 


10 9 


420 


-0 


.000 


054 


(218) 


AM 8LLF 


12 


-0.764 424 (213) 


1 


X 


10 9 


520 


-0 


.000 


255 


(465) 


AM 8LLG 


12 


-0.551 233 (148) 


1 


X 


10 9 


400 


0. 


.000 


265 


(686) 


AM 8tLfl 


6 


0.190 023 (111) 


1 


X 


10 9 


400 


+0 


.000 


747 


(691) 


AM 8U / 


12 


0.809 122 (217) 


1 


X 


10 9 


520 


-0 


.000 


484 


(529) 



for the gauge-invariant subgroups of group IV: 



i (8) 

IV(a) 

(8) 
2 IV(b) 

(8) 
2 IV(c) 

l (8) 
l IV(d) 



0.598 95 (14), 
0.822 49 (28), 
-1.138 91 (35), 
-0.990 72 (10). 



(51) 
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Summing up these terms one find that the contribution from all 180 diagrams of group IV 
is given by 

af^ = -0.708 19 (48). (52) 

VII. GROUP V DIAGRAMS 

There are 518 vertex diagrams without closed electron loops that contribute to Group 
V. We have treated them only in Version A thus far, in which these vertex diagrams are 
represented by 74 self-energy-like diagrams. Since some of them are time-reversal of others, 
the number of independent integrals is reduced to 47, which are shown in Figs. El and El 

Carrying out momentum space integration analytically by SCOONSCHIP and, more 
recently, by FORM, we obtain the magnetic moment projection of the Ward-Takahashi- 
summed amplitude as an integral over Feynman-parametric space of the form 

M 8a = — / (dz) orr9T .„ + orr o T/9 + 



256 J v ; \ 3U 2 V 3 3U 3 V 2 3U 4 V 
N + Z , N l + Z 1 , N 2 + Z 2 , N 3 + Z 3 



U 2 V 4 U 3 V 3 U 4 V 2 U 5 V J ' 1 ' 

where a = 1, ...,47. Eq through Z 3 are polynomials of scalar currents Bij and A{. See 
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301 ] for notation and definition. 
Standard subtractive renormalization relates the M 8a to the renormalized contribution 

by 

47 

ai v = a8c " ^ 

a=01 

where 

ag a = M$ a — renormalization terms. (55) 

Mg Q is the sum of rip vertex diagrams where np is given in the second column of Tables 
11X1 1X1 1X11 and 1X1 II For simplicity we denote these amplitudes as M a instead of M 8a in the 
following. To avoid possible confusion with M 2 , M 4 ,..., defined elsewhere, however, let us 
denote M x , M 2 , . . . , M 9 as M 01 , M 02 , . . . , M 09 . 

The two-step renormalization of eighth-order terms is in principle no more difficult than 
that of lower-order terms, except that it requires many more subtraction terms. UV and IR 

28 




M 22 M 23 M 24 



FIG. 8: Part of the electron self-energy-like diagrams Moi - M 2 4 representing (external- vertex- 
summed) integrals of subgroup V. 
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FIG. 9: Part of the electron self-energy-like diagrams M25 - M47 representing (external- vertex- 
summed) integrals of subgroup V. 
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subtraction terms of two-step renormalization scheme in Version A were originally derived 
manually. They were listed in Appendix D of Ref. (2^]. In view of the fact that the two-step 
renormalization scheme plays the central role in the success of our numerical integration 
approach, we have re-examined its validity by a FORM program. The FORM program is 



then translated into LATEX format 



Ezi- 



The result is listed in Appendix [B] together with 
he corresponding standard renormalization. We corrected some typographical errors in Ref. 
22|. 

Collecting contributions of all terms of Version A we obtain 

a { y ] = AM (8) - 5AM (6) A£ 2 

- AM (4) [4AL (4) + 3A5 (4) - 9(AB 2 ) 2 ] - AM^ASm^ 

- M 2 [2AL (6) + AE (6) - A<W 4) (4AL 2 * + AB 2 , - B 2 «[I})} 

- (M 2 » - M 2 »[/])[AcW 6) - AW 4) (5A5 2 + A5m 2 »)] 

+ M 2 AB 2 [10AL (4) + 6A5 (4) - 5(A5 2 ) 2 ], (56) 

where (see Ref. [22] for notation) 

47 H 
AM (8) = ^AM i5 AM (6) = ^AM 6a , 

i=01 a=A 

AAf( 4 *) = 2AM 4a(1 *) + AM 4a{2 *) + 2AM 4fc(1 *) + AM 4b(2 * h 

H 5 H 

AL(6) = ^$>aAL 6m , AB^ = J2vaAB 6a , 

a=A i=l a=A 
H 

A6m^ = J2v^6m 6a , (57) 

a=A 

where AMj represents the contribution of np vertices. AM 6a is also the sum of several 
vertex diagrams. Meanwhile, r] a in front of AL 6Qi , AB 6a , and A5mQ a is their multiplicity. 

Numerical evaluation of Group V integrals has been carried out by VEGAS. Tables HXl 
IXl IXI[ and IXHI list only the latest of long sequence of work carried out over more than 17 
years. Each integral was evaluated 5 to 10 times, starting with a relatively small number 
(typically ~ 10 7 ) of sampling points N$, increasing it gradually to more than 10 9 as more 
and more computing power became available. Although these early data are not included in 
the final statistics, they have played important roles in providing consistency check among 
different runs. In several cases two statistically independent runs of similar high precision 
have been combined statistically, as is indicated in the Tables. 
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It was seldom that VEGAS gave good result in the first try. In most cases we had to search 
for alternate approach and/or transform the integrand by various means, such as remapping, 
higher precision arithmetic, stretching, splitting, freezing, and their combinations. (See 
Appendix B of Ref. 3(| for terminology.) 

As is seen from Tables HX1 - IX11[ only ten integrals gave good results with real*8 algebra. 
The remaining 37 diagrams required real* 16 in some ways. In most cases splitting worked 
reasonably well. Integrals M 4, Mi 6 , Mi 8 , M 26 , M 28 , M 30 , M 42 , and M47 were evaluated in 
real* 16 without relying on splitting. M 4, M 16 , and M 18 were evaluated also in the adjustable 
precision 

Evaluation of M 40 in real* 16 ran into a rather puzzling problem: For relatively small 
values of N$ this integral seemed to behave reasonably well. As we increased N$, however, 
the error did not decrease as expected. For larger values of Ns, say Ns ~ 10 9 , for which each 
iteration now has small enough error bars, several independent sequences of runs gave results 
which are different from each other by much more than the error estimates generated by 
VEGAS. We suspected that this was caused by a particularly severe d-d problem which made 
it difficult to choose a reasonable mapping. To pinpoint the cause of this annoying problem 
we decided to run M 40 using real*32 arithmetic, which would reduce the d-d problem to a 
negligible level. Running M40 with Ns = 8 x 10 9 so that intrinsic resolution of individual 
iteration generated by VEGAS is ~ 0.003, we discovered that this integral fluctuates between 
two peaks which are far apart (~ 0.04) from each other. Once we discovered this problem, 
we went back to the beginning and chose a different mapping which does not have the 
double-peaking problem. Afterwards M40 behaved reasonably well. The value of M40 in 
Table IXTTl was obtained entirely in real*32. 

In comparison with M40 other diagrams were less troublesome in finding a reasonable 
mapping. However, diagram M 12 showed sign of considerable d-d problem. We therefore 
evaluated M i2 in real*32, too. The result is listed in Table HXl 

Combining the data from Tables HX| IX j IXl) and IXHI as well as auxiliary integrals defined 
in (|K7jl and listed in Table IXlTT] one can evaluate dy from (JHH|) . This leads to 

a p = -1.993 06 (343). (58) 

The precision we have reported here is more than adequate for comparison with experi- 
ment for the time being. However, to prepare for the confrontation with future measurements 
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of a e , further improvement of precision of some diagrams would be highly desirable. 
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APPENDIX A: ELIMINATION OF ALGEBRAIC ERROR OF GROUP V 



The number of independent integrals of Group V is 47 in Version A. They are denoted 



asMoi, M 02 ,---,M 47 



22j. All integrals of Group V diagrams can be constructed from just 



one "template" by permutation of tensor indices of photons. This implies that any error in 
one of the integrals will propagate to other integrals. On the other hand, verification of one 
integral will verify all other integrals of Group V. All integrands are created by the following 
procedures: 



33 



TABLE IX: Contributions of diagrams Mq\ — M13 of Fig. |HJ Some integrals are split into two parts: 
d-part is evaluated in real*8 and q-part is evaluated in real*16. a-part refers to the adjustable 
precision method developed by 48]. The superscript "a" means that the adjustable precision 
method is used in the whole domain. The superscript "q-d" on AM12 indicates that it is run 
with real*32. Integrals without such qualification are evaluated in double precision. Only the 
latest results are listed. AM04 is a combination of two independent computations with comparable 
statistical weights. 



Integral 


riF 


Value (Error) 


Sampling per 


No. of 






includin 


g n F 


iteration 


iterations 


AM 01 


7 


- 0.216 765 


(522) 






d-part 




-0.263 113 


(397) 


8 x 10 9 


139 


q-part 




0.046 348 


(338) 


2 x 10 8 


226 


AM 02 


14 


- 0.274 402 


(573) 






d-part 




- 0.223 509 


(347) 


8 x 10 9 


464 


q-part 




- 0.050 893 


(455) 


2 x 10 8 


320 


AM 03 


7 


- 0.681 633 


(518) 






d-part 




- 0.601 353 


(358) 


6 x 10 9 


332 


q-part 




- 0.080 280 


(374) 


2 x 10 8 


143 




14 


4.560 949 


(838) 


3.5 x 10 9 ,4 x 10 9 


233,209 


AM 05 


14 


2.397 618 


(266) 


4 x 10 9 


296 


AM 06 


14 


- 1.447 248 


(256) 


4 x 10 9 


286 


AM 07 


14 


0.065 257 


(282) 


4 x 10 9 


279 


AM 08 


14 


- 6.576 977 


(472) 






d-part 




- 6.724 530 


(343) 


1.6 x 10 10 


391 


q-part 




0.147 554 


(323) 


3 x 10 8 


125 


AM 09 


14 


0.937 503 


(447) 


8 x 10 9 


409 


AM10 


14 


17.859 588 


(600) 






d-part 




16.823 753 


(495) 


2 x 10 10 


403 


q-part 




1.035 835 


(338) 


4 x 10 8 


189 


AMn 


7 


2.433 695 


(336) 






d-part 




2.372 864 


(168) 


1 x 10 10 


378 


q-part 




0.060 831 


(291) 


1 x 10 8 


100 


AM[ q 2 - d) 


7 


- 4.096 653 


(871) 


1 x 10 9 


316 
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TABLE X: Contributions of M u - M 24 of Fig. El See the caption of Table HXI for notation. 
Only latest results are listed. AMis and q — parti of AMi6 are statistical combinations of two 
independent runs. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






includin 


g n F 


iteration 


iterations 


AM13 


14 


- 6.658 475 


(277) 


8 x 10 9 


280 


AMu 


14 


2.681 266 


(404) 






d-part 




2.653 346 


(289) 


7 x 10 9 


269 


q-part 




0.027 920 


(282) 


2 x 10 8 


114 


AM15 


14 


1.100 419 


(473) 






d-part 




1.199 683 


(298) 


5 x 10 9 


237 


q-part 




- 0.099 264 


(367) 


2 x 10 8 


236 


AMi 6 


14 


2.972 899 


(758) 






q — parti 




2.821 674 


(680) 


1 x 10 10 ,2 x 10 9 


150,799 


q - part 2 




0.151 225 


(334) 


4 x 10 8 


536 


AMn 


14 


3.157 998 


(415) 






d-part 




3.042 688 


(253) 


2.4 x 10 10 


189 


q-part 




0.115 310 


(329) 


1.2 x 10 8 


326 


AMfg ; 


14 


12.705 790 


(695) 


4 x 10 9 ,4 x 10 9 


188,145 


AM 19 


7 


- 0.834 502 


(266) 






d-part 




- 0.753 570 


(170) 


4 x 10 8 


221 


q-part 




- 0.080 932 


(204) 


2 x 10 7 


99 


AM 20 


14 


0.689 701 


(320) 






d-part 




0.604 560 


(299) 


1 x 10 9 


102 


q-part 




0.085 141 


(113) 


2 x 10 7 


105 


AM21 


7 


0.209 216 


(281) 






d-part 




0.191 592 


(245) 


4 x 10 8 


99 


q-part 




0.017 624 


(137) 


2 x 10 7 


100 


AM 22 


14 


- 0.281 374 


(454) 






d-part 




- 0.226 185 


(410) 


1 x 10 9 


110 


q-part 




- 0.055 189 


(194) 


2 x 10 7 


100 


AM 23 


14 


- 4.431 906 


(473) 






d-part 




-4.139 207 


(389) 


3 x 10 9 


289 


q-part 




- 0.292 699 


(268) 


4 x 10 7 


80 
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TABLE XI: Contributions of diagrams M 25 - M 39 of Fig. H See the caption of Table for 
notation. Only the latest results are listed. 



Integral 




Value (Error) 


Sampling per 


No. of 






lnrlnHin 


P~ Ti 


iteration 


iterations 


AMoa 


14 


2.206 446 


'404) 






d-part 




2.066 148 


'263) 


8 x 10 9 


209 


q-part 




0.140 276 


'314) 


2 x 10 8 


193 


AM 25 


7 


0.080 642 








a-part 




0.047 651 


lot) 


1 x 10 9 


100 


a-part 




0.032 991 


v o9oJ 


4 x 10 7 


106 


AM$ 


7 


1.829 535 


' V AV\ 

,545) 


9 x 10 8 , 1 x 10 9 


173,222 


AM 2 7 


14 


0.988 319 


^383) 






d-part 




n 000 A 1 T 

0.838 417 


v 25bJ 


4 x 10 9 


145 


a-part 




0.149 902 


v 2o5) 


1.2 x 10 s 


311 


AM^ q) 

Zo 


14 


/i nrn ^ frn 

— 4.950 479 


(690) 


2.6 x 10 9 ,4 x 10 9 


343,289 


AMoa 


7 


r\oo ann 

3.083 789 


(245) 






d-part 




3.023 232 


,130) 


6 x 10 9 


339 


a-part 




0.060 557 


(207) 


4 x 10 7 


234 




7 


— 3.721 347 


'630) 


2 x 10 9 , 2 x 10 9 


159 206 


AM,, 


7 


3.012 077 


179) 


2 x 10 9 


290 




14 


- 2.377 888 


(323) 






d-part 




i~\ Ail t \ A O 

— 2.343 943 


'284) 


4 x 10 8 


100 


G,-T)G,rt 

Vaj LJVaj 1 V 




- 0.033 945 


153) 


2 x 10 7 


100 




7 


- 1.194 130 


(250) 






nn.vt 

\AJ L/\AJ 1 V 




— 1.247 252 


'243) 


4 x 10 8 


98 


a-part 




0.053 122 


(57) 


2 x 10 7 


52 




14 


1.341 091 


(299) 






d-part 




1.237 527 


149) 


5 x 10 9 


411 


q-part 




0.103 564 


;259) 


4 x 10 7 


80 


AM35 


14 


- 0.539 732 


;384) 






d-part 




- 0.560 805 


;376) 


1 x 10 9 


100 


q-part 




0.021 073 


(75) 


2 x 10 7 


100 


AM 36 


14 


- 0.302 434 


(240) 






d-part 




- 0.252 129 


113) 


2 x 10 10 


275 


q-part 




- 0.050 305 


'211) 


4 x 10 7 


101 
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TABLE XII: Contributions of diagrams M 40 - M 47 of Fig. H See the caption of Table ED for 
notation. The superscript "q-d" on AM40 indicates that it is run with real*32. Only the latest 
results are listed. 



Integral 




Value (Error) 


Sampling per 


No. of 






m p n H 1 11 (T '>■) 7—1 


1 LCI tx LlUll 




AM37 


7 


0.487 701 (305) 


4 x 10 8 


90 


AM 38 


7 


— 3.094 936 (3ll) 


8 x 10 9 


147 


AM39 


14 


- 0.660 746 (290) 


4 x 10 9 


115 


A 71 /r{Q — d) 


14 


5.252 164 (977) 


1 x 10 9 


235 


AM41 


7 


- 1.416 949 (376) 






d-part 




- 1.107 727 (256) 


6 x 10 9 


129 


q-part 




- 0.309 222 (275) 


1.5 x 10 8 


245 


AM^ 


7 


- 4.859 600 (678) 


2 x 10 9 ,2 x 10 9 


168,226 


AM43 


7 


- 2.385 383 (221) 


1 x 10 9 


431 


AM44 


14 


3.468 870 (249) 






d-part 




3.288 775 (78) 


2 x 10 10 


441 


q-part 




0.180 091 (236) 


3 x 10 8 


232 


AM45 


7 


0.711 911 (411) 






d-part 




1.427 913 (217) 


4 x 10 9 


351 


q-part 




- 0.716 002 (349) 


3 x 10 8 


179 


AM 46 


7 


- 6.766 038 (242) 






d-part 




- 6.121 804 (86) 


3 x 10 10 


216 


q-part 




- 0.644 234 (226) 


6 x 10 7 


130 


am|? } 


7 


8.352 501 (585) 2.4 x 10 9 ,2 x 10 9 


205,179 



a. Momentum-space integration is carried out analytically by an algebraic manipulation 
program FORM |l5j . The result is expressed in terms of symbols A4 and (called scalar 
currents |22[ because they obey an analogue of Kirchhoff 's laws for electric current in which 
Feynman parameters play the role of resistance 0|). 

n 

b. Diagrams of Group V belong to five topologically distinct types |22|. Correspondingly, 
scalar currents Ai and Bij can be written explicitly as rational functions of Feynman param- 
eters using one of the five "templates" for Bij. Scalar currents must satisfy eight junction 
laws and four loop laws for each diagram. This imposes a very strict constraint on scalar 
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TABLE XIII: Auxiliary integrals for Group V. In this table we denote AM 4a * = 2AM 4a(1 „) + 
AM 4a ( 2 »). Similarly for AM 46 *. A LB d = Ya=i &L 6A i + \AB QA + 2A8m 6 A- Similarly for B LB d, 

etc. 

Integral Value (Error) Integral Value (Error) 



AB 2 * 


1.5 


AL 2 * 


-0.75 


B 2 *[I] 


-0.5 


A(5?T1-2* 


-0.75 


AM 4a * 


3.636 02 (11) 


AM 4fe « 


10.147 24 (14) 


AL 4 


0.465 024 (17) 


AB 4 


-0.437 094 (21) 


Albd 


1.370 25 (52) 


Blbd 


3.571 37 (26) 


Clbd 


-6.778 82 (39) 


2 x D LB d 


-8.783 15 (40) 


Elbd 


-1.592 21 (26) 


Flbd 


2.652 93 (36) 


2 x Glbd 


3.890 15 (33) 


Hlbd 


-0.232 26 (13) 


am( 6 ) 


4.471 013 (45) 







currents and provides powerful defense against trivial programming errors. 




FIG. 10: Chain diagrams obtained from eighth-order diagrams by omitting external lines and 
disregarding the distinction between electron and photon lines. The topological structure of chain 
diagram is crucial for determining the structure of "scalar currents" . 

c. Feynman-parametric integrals thus obtained are ultraviolet (UV) and/or infrared (IR) 
divergent in general. By a power counting rule designed for each divergence we can identify 
term by term which one is divergent [22]. Collection of these terms (in which only leading 
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terms are kept within each scalar current) forms a term required for renormalization. Its 
structure can be examined analytically by comparison with a known lower-order expression 
which is written as a Feynman-parametric integral over a function of its own scalar currents. 
These two integrals have different forms in general and their identity can be established only 
after some non-trivial analytic transformation. Comparison of these integrals provides an 
extremely powerful check of the subtraction term as well as major parts of mother integrand. 

d. The only remaining task is to show that scalar currents themselves have no error. This 
is not difficult since the number of scalar currents is not large. This is the primary reason 
for expressing the integrand in terms of scalar currents instead of Feynman parameters 
themselves. 

Initially, scalar currents were built by hands by two people working independently (21]. 
Comparison of resulting integrands enabled us to weed out human errors quickly. Recently 
we wrote programs in FORM (and MAPLE) which generate Bij for each chain diagrams of 
Fig. Eland show further that they satisfy all loop laws and junction laws. 

Around 1992 we re-derived the entire integrands by a computer program in a form dif- 
ferent from the original version. In this manner we obtained two sets of integrands with 
different algebraic structures. Their equivalence was checked by numerical spot check. This 
process was repeated for each renormalization term. Of course this comparison does not 
establish the validity of all terms. Since each diagram has several subtraction terms (as 
many as 20 in some cases), however, checks repeated for all subtraction terms cover a much 
larger fraction of the integrand than that of a single subtraction term and enhance strongly 
the confidence in the validity of the entire code. 

The procedure outlined above is especially potent for the diagrams that contain a second- 
order electron self-energy part. When the integrand is expressed as a function of scalar cur- 
rents, individual terms of the integrand, except those that contain the Feynman parameter 
z se for the electron of the self-energy part explicitly, are UV divergent. This means that 
the corresponding self-energy subtraction term, when expressed in terms of its own scalar 
currents, has a form identical with that of the mother integrand minus terms proportional to 
z se . Their difference comes solely from different structure of scalar currents for the mother 
and daughter integrands. 

This argument applies equally well to diagrams which contain two second-order electron 
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self-energy parts. In this case all terms of the mother integrand can be verified by applying 
the above argument to each self-energy part, verifying the correctness of the template com- 
pletely and thus all integrands obtained using the same template, namely, the entire Group 
V diagrams. In this manner we were able to confirm the validity of all terms of Group V 
integrals of Version A. Our very tight scheme a., b., c, and d leaves no room for further 
algebraic error. 



APPENDIX B: TWO-STEP RENORMALIZATION OF GROUP V INTEGRALS 



Following the discussion in Refs. 



121 . 1131 |22j, we adopted the two-step renormalization 



scheme in the numerical evaluation of g — 2. The K operation defined in Refs. [12, 22] picks up 
only the overall UV divergent part of the on-shell renormalization constant. It is used as the 
UV subtraction term in numerical calculation. The rest of the renormalization constants are 
combined in the end with those of other diagrams, which gives a finite contribution, referred 
to as the residual renormalization. The standard on-shell renormalization is the sum of K 
operation and residual renormalization. Only the renormalization scheme of Version A is 
given here since Group V has not been treated in Version B. 



1. Renormalization Scheme for Group V 

The eighth-order Group V diagrams have very complicated divergence structure. Only 
the K operation part was listed in Ref. because of page restriction. For convenience 
of access and completeness we list here K-renormalization and standard renormalization 
together. Few typographic errors in the previous version are corrected. 

We first list the standard on-shell renormalized amplitude dj expressed in terms of the 
unrenormalized integral Mj and various renormalizing terms and then present the interme- 
diate renormalization scheme which describes the relation between M$ and AM», the finite 
part of Mi defined by K operation. Substitution of Mi from the second part into the equa- 
tion for di expresses it in terms of AMj and pieces of renormalization constants as well as 
(unrenormalized) amplitudes of lower orders. Plugging them in (|54j) and using relations 
among lower order terms given in Appendix IB 21 we obtain Eq. (J56l) . 

M na is the unrenormalized amplitude of the nth-order diagram a, where n < 6. A quantity 
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with prefactor A implies that it is finite. The infrared subtraction terms are denoted as I na . 

Quantities L n , B n , and 5m n denote vertex renormalization constant, wave-function renor- 
malization constant, and mass renormalization constant of the nth-order, of the standard 
on-shell renormalization, respectively. In our approach by numerical integration, they are 
expressed as integrals over 3n, 3n — 1, and 3n — 1 Feynman parameters, which enable us 
to perform cancellation of divergence at every point in the domain of integration. In the 
parametric form the divergent terms of an integrand can be readily identified by means of a 
power counting rule for each divergence. The sum of these terms are denoted as L n , B n , and 
<5m n , and the remainder as L n , etc. For instance, L2 = L2 + L2, where L2 is UV-divergent 
part picked up by the K operation and L 2 is the residual UV-finite term. 

Within each subdiagram the electron lines are assumed to be consecutively labeled as [1, 
2, n m ], where n m = In for L n and In — 1 for B n and 5m n , and the line 1 is at the tip of 
the arrow attached to the electron line. 

The additional suffix such as c in L^ c indicates that it comes from the so-called corner 
diagram. Here we stick to the notation introduced previously. See Ref. |22j, p. 235, Fig. 8. 
The suffix (3*) means the insertion of two-point vertex (such as mass- vertex) in the fermion 
line 3. The primed quantity L^i') is the derivative amplitude obtained by applying — z\-^ 
operation on the integrand, where z\ is the Feynman parameter assigned to the fermion 
line 1. Note that L^nn is equal to L 4a; , but its UV-divergent part L± x ry\ is not equal to 
Li X . Since second order quantities such as B2 and Sm2 have only one electron line, it is not 
really necessary to distinguish different electron lines. We therefore use somewhat sloppy 
notations B2* and By instead of -62(1*) and -82(1') • Similarly, double insertions of two-point 
vertices are denoted as B 2 **, B 2 >*, B 2 ", etc. Since b~m,2 has no UV-finite part, 8ifi2 is replaced 
by 5rri2 everywhere. For L2, which contains electron lines 1 and 2, it is sometimes necessary 
to distinguish lines in which insertion is made. L 2 **x implies that two two-point vertices 
are inserted into the fermion line 1 of L2, while L 2 *t* means that one two-point vertex is 
inserted into the line 1 and another into the line 2. M^ a contains three electron lines 1,2,3 
and M 4a (!*.) means that two-point vertex insertion has been made twice in the electron line 
1, and so on. 

The validity of all equations has been confirmed algebraically by FORM. They are then 
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translated into LATEX format listed below 



47] 



On-Shell Renormalization for the Eighth-order Moments 

a i = Moi - 2M 6F L 2 - 2L 6F1 M 2 - 2L 4c M 4a + 6L 4c L 2 M 2 
+ 3M 4a (L 2 ) 2 -4(L 2 ) 3 M 2 

a 02 = M 02 - M 6D L 2 - M 6F B 2 - M &F{1 *)5m 2 - L 6D5 M 2 

+ L 4c ( 3 *)5m 2 M 2 - L 4c M 4b + L 4c 5m 2 M 2 * + 2L 4c B 2 M 2 - L 4s M 4a 

+ 2L 4s L 2 M 2 + M 4a{x *)bm 2 L 2 + M 4b (L 2 ) 2 + M 4a 5m 2 L 2 , + 2M 4a B 2 L 2 

- 25m 2 L 2 L 2 *M 2 - 5m 2 (L 2 ) 2 M 2 * - 3B 2 (L 2 ) 2 M 2 

a 03 = M 03 - 2M eD L 2 - M 6F B 2 - M GFm 5m 2 - 2L 6D1 M 2 

+ 2L 4c{1 * ) 5m 2 M 2 + 2L 4c B 2 M 2 + 2L 4s L 2 M 2 + 2M 4a(l * ) 5m 2 L 2 + M 4b (L 2 ) 2 
+ 2M 4a B 2 L 2 - 25m 2 L 2 L 2 ,M 2 - 5m 2 (L 2 ) 2 M 2 « - 3B 2 (L 2 ) 2 M 2 

a 04 = M 04 - M 6A L 2 - 2M 6D B 2 - M 6D (i*)drn 2 - M 6D ( 3 *)5m 2 

- L 6A1 M 2 + 2L 4s{lf) 5m 2 M 2 + 2L 4s B 2 M 2 + 2M 4a(1 * ) 6m 2 B 2 + 2M 4b{1 , ) 5m 2 L 2 
+ 2M 4b B 2 L 2 + M 4a {B 2 ) 2 + M 4a(in * } (5m 2 ) 2 - 25m 2 B 2 L 2 M 2 , - 25m 2 B 2 L 2 ,M 2 

- (5m 2 ) 2 L 2 M 2 „ - (5m 2 ) 2 L 2 * n M 2 - 2{B 2 ) 2 L 2 M 2 

a 5 = M 05 - M 6H L 2 - L 6H1 M 2 - L 6F2 M 2 - L 4x M 4a 
+ 3L 4x L 2 M 2 

o-oe = M oe — M 6F L 2 — M 6G L 2 — L 6F3 M 2 — L 6G5 M 2 

+ 3L 4c L 2 M 2 - L 4l M 4a + 2L 4l L 2 M 2 + 2M 4a (L 2 ) 2 - 3(L 2 ) 3 M 2 

007 = M 07 — M SF L 2 — MqqL 2 — Lqg\M 2 — L§ F)2 M 2 

- L 4c M 4a + 4L 4c L 2 M 2 + L 4l L 2 M 2 + 2M 4a (L 2 ) 2 - 3(L 2 ) 3 M 2 



42 



a s = M 08 - Mq C L 2 - 2M 6D L 2 - L 6C iM 2 - 5m 4a M 4a(1 ») 

+ 5m 4a L 2 M 2 * + 5m 4a L 2 *M 2 - B 4a M 4a + 2B 4a L 2 M 2 + 2L 4s L 2 M 2 

+ 2M 4a(1 * ) 5m 2 L 2 + 2M 4b (L 2 ) 2 + 2M 4a B 2 L 2 - 25m 2 L 2 L 2 *M 2 - 25m 2 (L 2 f M 2 , 

- AB 2 {L 2 fM 2 

a g = M og - M 6E L 2 - M 6F B 2 - M &F ( 2 *)5m 2 - L 6E iM 2 

- L 6D3 M 2 + L 4ci i* ) 8m 2 M 2 + L 4c{2 * ) 5m 2 M 2 + 2L 4c B 2 M 2 - L 4s M 4a 

+ 3L 4s L 2 M 2 + M 4a{2 ^m 2 L 2 + M 4a 5m 2 L 2 * + 2M 4a B 2 L 2 - 35m 2 L 2 L 2 *M 2 

- 3B 2 {L 2 ) 2 M 2 

a w = Mio - M 6B L 2 - M 6D B 2 - M %D{2 *)5m 2 - L 6B1 M 2 

+ L 4s{2 *)5m 2 M 2 - 5m 4h M 4a{ i*) + 5m 4b L 2 M 2 * + 5m 4b L 2 *M 2 - B 4b M 4a 

+ 2B 4b L 2 M 2 + L 4s B 2 M 2 + M 4a (i*)<5m 2 <5m 2 * + M 4a ( 1 *- ) 5m 2 B 2 + M 4b ( 2 *- ) 5m 2 L 2 

+ M 4b B 2 L 2 + M 4a 5m 2 B 2 * + M 4a {B 2 ) 2 — 5m 2 5m 2 * L 2 M 2 * — Sm 2 Sm 2 *L 2 *M 2 

- 5m 2 B 2 L 2 M 2 * - bm 2 B 2 L 2 ,M 2 - 25m 2 B 2 *L 2 M 2 - 2(B 2 ) 2 L 2 M 2 

a n = Mn - 2M 6D B 2 - 2M 6Dm 5m 2 - 2L 4s M 4b + 2L 4s 5m 2 M 2 * 

+ 2L 4s B 2 M 2 + 2M 4a (i*)5m 2J B2 + 2M 4b 5m 2 L 2 * + 2M 4b B 2 L 2 + M 4a (B 2 ) 2 

+ M 4a(1 ,3*)(5m 2 ) 2 - 25m 2 B 2 L 2 M 2 * - 25m 2 B 2 L 2 *M 2 - 2{bm 2 ) 2 L 2 »M 2 » - 2{B 2 ) 2 L 2 M 2 

a 12 = M 12 - 3M %A B 2 - 2M &A{1 *)5m 2 - M 6j4 ( 3 *)5m 2 + QM 4b ^*- ) 5m 2 B 2 

+ 3M 46 (5 2 ) 2 + 2M 46(in ,)(5m 2 ) 2 + M 46(1 , 3 »)((5m 2 ) 2 - 35m 2 (B 2 ) 2 M 2 * - 3(5m 2 ) 2 B 2 M 2 ** 

- (5m 2 ) 3 M 2 *** - (5 2 ) 3 M 2 

ai 3 = Mi 3 - M m B 2 - M m{1 *)5m 2 - L 6D4 M 2 + L 4x{4 *)bm 2 M 2 

- L 4x M 4b + L 4 Jm 2 M 2 * + 2L 4x B 2 M 2 

a 14 = M M - M 6D L 2 - M 6G B 2 - M 6G(5 *)5m 2 - L 6D3 M 2 

+ L 4c(1 *)5m 2 M 2 + L 4c B 2 M 2 + L 4s L 2 M 2 - L 4 iM 4b + L 4l 5m 2 M 2 * 

+ L 4l B 2 M 2 + M 4a{1 * } 5m 2 L 2 + M 4b {L 2 f + M 4a B 2 L 2 - 5m 2 L 2 L 2 *M 2 

- 5m 2 (L 2 ) 2 M 2 , - 2B 2 [L 2 ) 2 M 2 
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a 15 = Mi 5 - Mq D L 2 - M 6G B 2 - M &G( i*)8m 2 - L 6A2 M 2 

+ L 4Kl * ) 5m 2 M 2 - L 4c M 4b + L 4c 8m 2 M 2 , + L 4c B 2 M 2 + L 4s L 2 M 2 
+ L 4l B 2 M 2 + M 4a{1 * } Sm 2 L 2 + M 4b (L 2 f + M 4a B 2 L 2 - Sm 2 L 2 L 2 *M 2 

- Sm 2 (L 2 ) 2 M 2 * - 2B 2 (L 2 ) 2 M 2 

ai6 = Mi 6 - 2M 6 a-^2 - M 6C B 2 - M 6 c(i*)Sm 2 - 8m 4a M 4b{1 *) 

+ <5m 4a 5m 2 M 2 ** + 8m 4a B 2 M 2 * - B 4a M 4b + B 4a 5m 2 M 2 * + B 4a B 2 M 2 

+ 4:M 4b{1 *)5m 2 L 2 + AM 4b B 2 L 2 - 48m 2 B 2 L 2 M 2 * - 2{8m 2 ) 2 L 2 M 2 „ - 2{B 2 ) 2 L 2 M 2 

an = M 17 - M 6D B 2 - M 6E B 2 - M 6E ( 1 *- ) 5m 2 - M 6D(4 »)5m 2 

- L 6A3 M 2 + 2L 4s{1 * ) 5m 2 M 2 - L 4s M 4b + L 4s 5m 2 M 2 * + 3L 4s B 2 M 2 

+ M M1 * } 5m 2 B 2 + M 4a{2 * } 5m 2 B 2 + M 4b 8m 2 L 2 * + M 4b B 2 L 2 + M 4a (B 2 f 

+ M 4a (i* 2 *)(5m 2 ) 2 - 5m 2 B 2 L 2 M 2 * - 35m 2 B 2 L 2 *M 2 - (5m 2 ) 2 L 2 *M 2 * - (8m 2 ) 2 L 2 *f*M 2 

- 2(B 2 ) 2 L 2 M 2 

a 18 = Mis - M 6A B 2 - M 6B B 2 - M (jB{l , ) 8m 2 - M &A{2 * ) 8m 2 

- 5m 4b M 4b ^) + 5m 4b 5m 2 M 2 ** + 5m 4b B 2 M 2 * - B 4b M 4b + B 4b 8m 2 M 2 * 

+ B 4b B 2 M 2 + M ib ^8m 2 8m 2 * + 2M 4b ( 1 *- ) 8m 2 B 2 + M 4b ( 2 *)8m 2 B 2 + M 4b 8m 2 B 2 * 

+ 2M 4b {B 2 ) 2 + M 46(1 * 2 *)(5m 2 ) 2 - 8m 2 8m 2 *B 2 M 2 * - Sm 2 B 2 B 2 *M 2 - 28m 2 (B 2 ) 2 M 2 * 

- (8m 2 ) 2 8m 2 *M 2 ** - (8m 2 ) 2 B 2 M 2 ** - (8m 2 ) 2 B 2 *M 2 * - (B 2 ) 3 M 2 

dig = Mi 9 - 2L 6H2 M 2 

a 20 = M 20 - MqhL 2 - L 6F2 M 2 - L 6G4 M 2 + 2L 4x L 2 M 2 
a 21 = M 21 - 2L 6G2 M 2 

a 22 = M 22 - Mq G L 2 - L 6Fl M 2 - L %C2 M 2 - L 4c M 4a 
+ 3L 4c L 2 M 2 + L 4i L 2 M 2 + M 4a (L 2 ) 2 - 2(L 2 ) 3 M 2 
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a 23 = M 23 - M m B 2 - M mi 2*)5m 2 - L 6E2 M 2 - L 6D4 M 2 
+ L 4x{ i,)5m 2 M 2 + L 4x{2 *)5m 2 M 2 + 2L 4x B 2 M 2 

a 24 = M 24 - M 6G B 2 - M 6G ( 2 *)5m 2 - L 6B2 M 2 - L 6D5 M 2 

+ L 4c{ i*)5m 2 M 2 + L Al{2 *)8m 2 M 2 + L 4c B 2 M 2 - L 4s M 4a + 2L 4s L 2 M 2 
+ L 4l B 2 M 2 + M 4a bm 2 L 2 » + M 4a B 2 L 2 - 25m 2 L 2 L 2 *M 2 - 2B 2 {L 2 ) 2 M 2 

a 25 = M 25 - 2M 6G L 2 - 2L 6D2 M 2 + 2L 4c L 2 M 2 + 2L 4l L 2 M 2 
+ M 4a (L 2 ) 2 -2(L 2 ) 3 M 2 

a 26 = M 26 - 2M 6C L 2 - 5m 6F M 2 * - B 6F M 2 + 25m 4a L 2 M 2 * 

+ 2B 4a L 2 M 2 - 2L 4c M 4b + 2L 4c 5m 2 M 2 * + 2L 4c B 2 M 2 + 3M 4b (L 2 ) 2 

- 35m 2 {L 2 ) 2 M 2 * - 3B 2 [L 2 ) 2 M 2 

a 27 = M 27 - M 6E L 2 - M 6G B 2 - M 6G(4 .)5m 2 - L 6D1 M 2 

- L 6A2 M 2 + L 4c{1 *)5m 2 M 2 + L 4U i*)5m 2 M 2 + L 4c B 2 M 2 + 2L 4s L 2 M 2 

+ L 4l B 2 M 2 + M 4a(2 *)5m 2 L 2 + M 4a B 2 L 2 - 25m 2 L 2 L 2 *M 2 - 2B 2 (L 2 ) 2 M 2 

a 28 = M 28 - M 6B L 2 - M 6G B 2 - M &G{2 *)5m 2 - 5m 6D M 2 * 

- B 6D M 2 + 8m 4a (i*)8m 2 M 2 * + B 4a ^*)5m 2 M 2 + 5m 4a B 2 M 2 * + 5m 4b L 2 M 2 * 
+ B 4a B 2 M 2 + B 4b L 2 M 2 - L 4s M 4b + L 4s 5m 2 M 2 * + L 4s B 2 M 2 

+ M 4b ( 2 .)5m 2 L 2 + M 4b 5m 2 L 2 * + 2M 4b B 2 L 2 — 5m 2 5m 2 * L 2 M 2 * — 25m 2 B 2 L 2 M 2 * 

- 5m 2 B 2 L 2 *M 2 - 5m 2 B 2 *L 2 M 2 - (5m 2 ) 2 L 2 , M 2 * - 2{B 2 ) 2 L 2 M 2 

a 29 = M 29 - 2M eE B 2 - 2M 6E{2 *)5m 2 - 2L 6M M 2 + AL 4s{l *)5m 2 M 2 

+ 4L 4s B 2 M 2 + 2M 4a{2 * ) 5m 2 B 2 + M 4a [B 2 ) 2 + M 4a(2 , 2 *)(5m 2 ) 2 - 45m 2 B 2 L 2 *M 2 

- 2(5m 2 ) 2 L 2 ** t M 2 - 2{B 2 fL 2 M 2 

a 30 = M 30 - 2M 6B B 2 - 2M &B{2 ,)5m 2 - 5m &A M 2 , - B 6A M 2 

+ 25m 4b{1 * ) 5m 2 M 2 * + 2 J B 46(1 * ) <5m 2 M 2 + 25m 4b B 2 M 2 * + 2B 4b B 2 M 2 + 2M 4b{2 , ) 5m 2 B 2 
+ M 4b (B 2 ) 2 + M 4fe(2 * 2 .)(5m 2 ) 2 - 25m 2 8m 2 *B 2 M 2 * - 25m 2 B 2 B 2 *M 2 - 5m 2 (B 2 ) 2 M 2 * 

- (5m 2 ) 2 5m 2 **M 2 * - (5m 2 ) 2 B 2 **M 2 - (B 2 ) 3 M 2 
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a 31 = M 3 i - 2L 6H3 M 2 

a 32 — M32 — L§h 2 M 2 — L eG3 M 2 

a 33 = M 33 - 2L 6G3 M 2 

a 34 = M 34 - L mi M 2 - L 6C3 M 2 - L 4x M 4a + 2L 4x L 2 M 2 
035 = M 35 - M m L 2 - L 6E3 M 2 - L 6G4 M 2 + 2L 4x L 2 M 2 

a 36 = M 36 - M 6G L 2 - L 6B3 M 2 - L 6G5 M 2 + L 4c L 2 M 2 

- L 4l M 4a + 3L 4l L 2 M 2 + M 4a (L 2 ) 2 - 2(L 2 ) 3 M 2 

a 37 = M 37 - 2L 6G2 M 2 

a 38 = M 38 - 5m 6H M 2 * - 5 6H M 2 - 2L 4x M 4b + 2L 4x 5m 2 M 2 * + 2L 4x B 2 M 2 

a 39 = M 39 - M 6G L 2 - L 6G1 M 2 - L 6C2 M 2 - L 4c M 4a 
+ 3L 4c L 2 M 2 + L 4l L 2 M 2 + M 4a {L 2 f - 2(L 2 ) 3 M 2 

a 40 = M 40 - M ec L 2 - 5m 6G M 2 * - B 6G M 2 + 5m 4a L 2 M 2 * 
+ B 4a L 2 M 2 - L 4c M 4b + L 4c 5m 2 M 2 * + L 4c B 2 M 2 - L 4i M 4b 
+ L 4l 5m 2 M 2 * + L 4l B 2 M 2 + 2M 4b (L 2 ) 2 - 25m 2 (L 2 ) 2 M 2 * - 2B 2 {L 2 ) 2 M 2 

a 41 = M 4 i - 2M 6E L 2 - 2L eG1 M 2 - Sm 4a M 4a{2 , } + 25m 4a L 2 *M 2 

- B 4a M 4a + 2B 4a L 2 M 2 + AL 4s L 2 M 2 + 2M4 a(2 ,)5m 2 L 2 + 2M 4a B 2 L 2 

- A5m 2 L 2 L 2 ,M 2 - AB 2 (L 2 ) 2 M 2 

a 42 = M 42 - 2M 6B L 2 - <5m 6C M 2 * - B 6C M 2 - 5m 4a M 4b{2 ,) 

+ 5m 4a 5m 2 * M 2 * + 5m 4a B 2 *M 2 + 25m 4b L 2 M 2 * - B 4a M 4b + B 4a Sm 2 M 2 * 
+ B 4a B 2 M 2 + 2B 4b L 2 M 2 + 2M 4f >( 2 *)5m 2 L 2 + 2M 4b B 2 L 2 - 25m 2 <5m 2 *L 2 M 2 * 

- 25m 2 B 2 L 2 M 2 * - 25m 2 B 2 *L 2 M 2 - 2(B 2 ) 2 L 2 M 2 



46 



a 43 = M 43 - M m B 2 - M 6H ( 3 *)5m 2 - 2L 6E2 M 2 + 2L 4x{2 * ) 5m 2 M 2 
+ 2L 4x B 2 M 2 



a 44 = M 44 - M 6G B 2 - M 6G ( 3 *)5m 2 - L 6E1 M 2 - L 6B2 M 2 

+ L Ac ( 2 *)5m 2 M 2 + L Al ( 2 *- ) 5m 2 M 2 + L 4c B 2 M 2 - L 4s M 4a + 2L 4s L 2 M 2 
+ L 4l B 2 M 2 + M 4a 5m 2 L 2 , + M 4a B 2 L 2 - 25m 2 L 2 L 2 *M 2 - 2B 2 (L 2 ) 2 M 2 

a 45 = M 45 - M 6C B 2 - M 6C (3*)Sm 2 - 5m 6E M 2 * - B eE M 2 

+ 5m 4a ( 2 *- ) 5m 2 M 2 * + B 4a ( 2 *- ) 5m 2 M 2 + 5m 4a B 2 M 2 * + B 4a B 2 M 2 - 2L 4s M 4b 
+ 2L 4s 5m 2 M 2 * + 2L 4s B 2 M 2 + 2M 4b bm 2 L 2 » + 2M 4b B 2 L 2 - 2bm 2 B 2 L 2 M 2 » 

- 25m 2 B 2 L 2 *M 2 - 2(5m 2 ) 2 L 2 * M 2 * - 2(B 2 ) 2 L 2 M 2 

a 46 = M 46 - Mq E B 2 - M 6E ( 3 *)5m 2 - 2L 6B1 M 2 + 2L 4s{2 ,)5m 2 M 2 

- 5m 4& M 4a (2*) + 25m 4b L 2 *M 2 - B 4b M 4a + 2B 4b L 2 M 2 + 2L 4s B 2 M 2 

+ M 4a{2 *)5m 2 5m 2 * + M 4a{2 *)5m 2 B 2 + M 4a 5m 2 B 2 * + M 4a (B 2 ) 2 - 25m 2 5m 2 * L 2 * M 2 

- 2bm 2 B 2 L 2 ,M 2 - 2bm 2 B 2 *L 2 M 2 - 2(B 2 ) 2 L 2 M 2 

a 47 = M 47 - M 6B B 2 - M eB (3*)8m 2 - 8m eB M 2 * - B 6B M 2 

+ 5m 4b ( 2 ,)5m 2 M 2 * + B 4b{2 ,)5m 2 M 2 - 5m 4b M 4b{2 *) + 5m 4b 5m 2 *M 2 * + 5m 4b B 2 M 2 * 

+ bm 4b B 2 »M 2 - B 4b M 4b + B 4b bm 2 M 2 » + 2B 4b B 2 M 2 + M 4b{2 *)5m 2 5m 2 * 

+ M 4b (2*)Sm 2 B 2 + M 4b 5m 2 B 2 * + M 4b (B 2 ) 2 — 5m 2 8m 2 * B 2 M 2 * — 5m 2 5m 2 * B 2 * M 2 

- 5m 2 (5m 2 *) 2 M 2 * - 25m 2 B 2 B 2 *M 2 - 5m 2 (B 2 ) 2 M 2 * - (5m 2 ) 2 B 2 *M 2 * - (B 2 ) 3 M 2 

Intermidiate Renormalization for the Eighth-order Moments 

M i = AM 01 + 2M 6F L 2 + 2L 6F1 M 2 + 2L 4c M 4a - QL 4c L 2 M 2 
-3M 4a (L 2 ) 2 + A(L 2 ) 3 M 2 
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M 02 = AM 02 + M 6D L 2 + M QF B 2 + M 6F{1 * ) 5m 2 + h F \M 2 
+ L 6D5 M 2 + L 4c M 4b - L 4c 5m 2 M 2 * - L 4c B 2 M 2 + L 4s M 4a 

- 2L 4s L 2 M 2 - L 4c{v) B 2 M 2 + AL 4c M 2 I 2 - M Ml * ) 5m 2 L 2 - M 4b (L 2 ) 2 

- M 4a B 2 L 2 - M 4a B 2 L 2 , + 5m 2 {L 2 fM 2 * + 2B 2 L 2 M 2 L 2 , + B 2 (L 2 ) 2 M 2 

M 03 = AM 03 + 2M 6D L 2 + M 6F B 2 + M 6F(3 * ) 5m 2 + 2L 6D1 M 2 

+ I 6F3 M 2 - 2L 4s L 2 M 2 - 2L 4c({VY) B 2 M 2 - 2M 4a{1 * ) 5m 2 L 2 - M 4b {L 2 f 

- 2M 4a B 2 L 2 + 5m 2 (L 2 ) 2 M 2 * + 2B 2 L 2 M 2 L 2 , + B 2 (L 2 ) 2 M 2 

M 04 = AM 04 + M 6A L 2 + 2M 6D B 2 + M 6D{1 * ) 5m 2 + M 6D{3 « ) 5m 2 
+ L 6A1 M 2 + I 6D1 M 2 + I 6D3 M 2 — L 4s ^ v y ) B 2 M 2 — L 4s ( 3 ,- ) B 2 M 2 

- 2M 4a{1 * } 5m 2 B 2 - 2M 4bil * } 5m 2 L 2 - 2M 4b B 2 L 2 - M 4a {B 2 ) 2 - M 4a(in *)(5m 2 ) 2 
+ L 2 ,,{B 2 ) 2 M 2 + 2bm 2 B 2 L 2 M 2 , + {5m 2 ) 2 L 2 M 2 ** + {B 2 ) 2 L 2 M 2 

M 05 = AM 05 + M m L 2 + L mi M 2 + L 6F2 M 2 + L 4x M 4a 

- 3L 4x L 2 M 2 

M 06 = AM 06 + M 6F L 2 + M 6G L 2 + L 6F3 M 2 + L 6G5 M 2 

- 3L 4c L 2 M 2 + L 4l M 4a - 2L 4l L 2 M 2 - 2M 4a (L 2 ) 2 + 3(L 2 ) 3 M 2 

M 07 = AM 07 + M 6F L 2 + M 6G L 2 + L 6G1 M 2 + L 6D2 M 2 

+ L 4c M 4a - AL 4c L 2 M 2 - L 4l L 2 M 2 - 2M 4a (L 2 ) 2 + 3(L 2 ) 3 M 2 

M 08 = AM 08 + M 6C L 2 + 2M 6D L 2 + L 6C1 M 2 + 5m 4a M 4a{1 * } 

- 5m 4a L 2 M 2 * + B 4a M 4a - B 4a L 2 M 2 - B 4a M 2 L 2 , - 2L 4s L 2 M 2 

+ A5m 4a / 4a( i.) - 2M 4a(l * ) 5m 2 L 2 - 2M 4b [L 2 ) 2 - 2M 4a B 2 L 2 + AM 4a / 4c 
+ 25m 2 (L 2 ) 2 M 2 * + 2B 2 L 2 M 2 L 2 , + 2B 2 {L 2 ) 2 M 2 

M 09 = AM 09 + M 6E L 2 + Mq F B 2 + M 6F{2 * ) 5m 2 + L 6E1 M 2 

+ / 6 f 2 M 2 + L GD3 M 2 + L 4s M 4a - 3L 4s L 2 M 2 - L 4c(v) B 2 M 2 

- L 4c{2l) B 2 M 2 - M 4a{2 * } 5m 2 L 2 - M 4a B 2 L 2 - M 4a B 2 L 2 , + 3B 2 L 2 M 2 L 2 , 
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Mio = AMio + M 6B L 2 + M 6D B 2 + M 6D{2 * ) 5m 2 + L QB1 M 2 

+ h D 2M 2 + 5m 46 M 4a(1 *) - 5m 4b L 2 M 2 * + B 4b M 4a - B 4b L 2 M 2 

- B 4b M 2 L 2 > + A5m 46 / 4a (i*) - L is{2 ')B 2 M 2 - M 4a ^*)5m 2 5m 2 * - M 4a ^*)B 2 5m 2 > 

- M 4b(2 * } 5m 2 L 2 - M 4b B 2 L 2 - M 4a B 2 B 2 , + AM 4b I 4c + h c M 2 I 2 
+ 5m 2 L 2 8m 2 * M 2 * + B 2 L 2 M 2 *5m 2 i + B 2 L 2 M 2 B 2 i + B 2 M 2 B 2 iL 2 i 

M n = AM n + 2M 6D B 2 + 2M 6D(5 * ) 5m 2 + 2I 6m M 2 + 2L 4s M 4b 

- 2L 4s 5m 2 M 2 * - 2L 4s B 2 M 2 + 2AL 4s M 2 I 2 - 2M 4a(1 * ) 5m 2 B 2 - 2M 4b B 2 L 2 , 

- M 4a {B 2 f - M 4a(1 * 3 ,)(5m 2 ) 2 + 28m 2 B 2 M 2 *L 2 , + 2{B 2 ) 2 M 2 L 2 , 

M12 = AMi 2 + 3M 6A B 2 + 2M 6A(1 »)5m 2 + M 6A( 3*)5m 2 + 2I 6Al M 2 

+ I 6A3 M 2 - 6M 4b{in Sm 2 B 2 - 3M 4b {B 2 f - 2M 4b(1 * 1 * ) (5m 2 ) 2 - M^^Sm^ 2 
+ 35m 2 (B 2 ) 2 M 2 * + 3(5m 2 ) 2 B 2 M 2 ** + (5m 2 ) 3 M 2 *** + (B 2 ) 3 M 2 

Mi 3 = AM 13 + M 6//J B 2 + M m{1 * ) 5m 2 + I mi M 2 + L 6D4 M 2 

+ L 4x M 4b — L 4x 5m 2 M 2 * — L 4x B 2 M 2 + AL 4x M 2 I 2 — L 4x (i^ B 2 M 2 

M 14 = AM 14 + MqdL 2 + M 6G B 2 + M 6G( 5,)<5m 2 + L 6D3 M 2 
+ h Gb M 2 - L 4s L 2 M 2 + L 4[ M 4b - L 4l 5m 2 M 2 * - L 4i B 2 M 2 
+ AL 4l M 2 I 2 - L 4c(v) B 2 M 2 - M 4a(l *)5m 2 L 2 - M 4b [L 2 ) 2 - M 4a B 2 L 2 
+ 5m 2 (L 2 ) 2 M 2 * + B 2 L 2 M 2 L 2 , + B 2 (L 2 ) 2 M 2 

M 15 = AM 15 + M 6Z) L 2 + M 6Gj B 2 + M 6G(1 »)<5m 2 + J 6G iM 2 

+ L 6 A2M 2 + L 4c M 46 - L 4c 5m 2 M 2 * - L 4c B 2 M 2 - L 4s L 2 M 2 

+ AL 4c M 2 I 2 - L 4l(v) B 2 M 2 - M 4a(l «)5m 2 L 2 - M 4b [L 2 ) 2 - M 4a B 2 L 2 

+ 5m 2 (L 2 ) 2 M 2 * + B 2 L 2 M 2 L 2 , + B 2 (L 2 ) 2 M 2 

M 16 = AM 16 + 2M 6j4 L 2 + M 6C B 2 + M 6C(1 * ) 5m 2 + / 6C1 M 2 

+ 1/2J 6C M 2 + 5m 4a M 4f ,(i») - 5m 4a 5m 2 M 2 ** - 5m 4a B 2 M 2 * + B 4a M 4b 

- B 4a 5m 2 M 2 * - B 4a B 2 M 2 + A8m 4a I 4b (i*) - A5m 4a M 2 I 2 * - 4M 4b ^)5m 2 L 2 

- AM 4b B 2 L 2 + AM 4a I 4s + A5m 2 B 2 L 2 M 2 * + 2{bm 2 ) 2 L 2 M 2 »» + 2{B 2 ) 2 L 2 M 2 
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M 17 = AM 17 + M 6D B 2 + M QE B 2 + M 6E(1 * ) 5m 2 + M 6D(4 *)5m 2 
+ I 6E1 M 2 + L 6A3 M 2 + I 6D4 M 2 + L 4s M 4b - L 4s 5m 2 M 2 * 

- L 4s B 2 M 2 + AL 4s M 2 I 2 - 2L 4s(v) B 2 M 2 - M 4a{1 * } 5m 2 B 2 - M 4a(2 * )( 5m 2J B 2 

- M 4h B 2 L 2 , - M 4a {B 2 f - M Ml * 2 * } (5m 2 ) 2 + L 2 n(B 2 ) 2 M 2 + 5m 2 B 2 M 2 *L 2 , 
+ (B 2 ) 2 M 2 L 2 , 

Mis = AM 18 + M 6Aj B 2 + M 6Bj B 2 + M 6B{1 * ) 5m 2 + M 6A(2 * } 5m 2 

+ UbiM 2 + Ua2M 2 + l/2J 6jB M 2 + 8m 4b M 4b (i*) - 5m 4b 5m 2 M 2 ** 

- 5rh 4b B 2 M 2 * + B 4b M 4b - B 4b 5m 2 M 2 * - B 4b B 2 M 2 + ASm 4b I 4b ^) 

- A5m 4b M 2 I 2 * - M 4b{ i*)5m 2 B 2 - M 4b{1 *)5m 2 5m 2 * - M 4b ^*)B 2 5m 2 t - M 4b{2 *)5m 2 B 2 

- M 4b B 2 B 2 > - M 4b {B 2 f - M 4fe( i* 2 *)(5m 2 ) 2 + AM 4b I 4s + I 4s M 2 I 2 

+ 5m 2 B 2 8m 2 * M 2 * + 5m 2 B 2 M 2 *B 2 i + 5m 2 B 2 M 2 **5rh 2 t + (5m 2 ) 2 5rh 2 *M 2 ** + (B 2 ) 2 M 2 *5rh 2 / 
+ [B 2 ) 2 M 2 B 2 , 

M19 = AM 19 + 2L m2 M 2 

M 20 = AM 20 + M m L 2 + L 6F2 M 2 + L 6G4 M 2 - 2L 4x L 2 M 2 
M 21 = AM 2l + 2L 6G2 M 2 

M 22 = AM 22 + M 6G L 2 + L 6F1 M 2 + L 6G2 M 2 + L 4c M 4a 

- 3L 4c L 2 M 2 - L 4l L 2 M 2 - M 4a {l 2 f + 2(L 2 ) 3 M 2 

M 23 = AM 23 + M m B 2 + M &H{2 , ) 5m 2 + L 6E2 M 2 + I m2 M 2 
+ L 6D4 M 2 — L 4x (i>)B 2 M 2 — L 4x ( 2 ')B 2 M 2 

M 24 = AM 24 + M 6G B 2 + M 6G(2 * ) 5m 2 + L &B2 M 2 + I 6G2 M 2 

+ L 6D5 M 2 + L 4s M 4a - 2L 4s L 2 M 2 - L 4c(y) B 2 M 2 - L 4l{2>) B 2 M 2 

- M 4a B 2 L 2 , + 2B 2 L 2 M 2 L 2 , 
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M 25 = AM 25 + 2M 6G L 2 + 2L eD2 M 2 - 2L 4c L 2 M 2 - 2L 4l L 2 M 2 
-M 4a (L 2 ) 2 + 2(L 2 ) 3 M 2 

M 26 = AM 26 + AM 6F I 2 + 2M 6C L 2 + 5m 6F M 2 * + B 6F M 2 

+ A5m 6F M 2 *[I} - 25m 4a L 2 M 2 * - 2B 4a L 2 M 2 + 2L 4c M 4b - 2L 4c 5m 2 M 2 * 

- 2L 4c B 2 M 2 - 3M 4b {L 2 ) 2 + 35m 2 {L 2 ) 2 M 2 * + 3B 2 [L 2 ) 2 M 2 

M 27 = AM 27 + M 6E L 2 + M 6G B 2 + M 6G{4 * ) 5m 2 + L 6D1 M 2 

+ Lqa 2 M 2 + Iq G4 M 2 — 2L 4s L 2 M 2 — L 4c ^yy)B 2 M 2 — L 4 \( V )B 2 M 2 

- M 4a{2 *)5m 2 L 2 - M 4a B 2 L 2 + 2B 2 L 2 M 2 L 2 , 

M 2S = AM 28 + AM 6D I 2 + M 6B L 2 + M 6C B 2 + M 6C{2 * ) 5m 2 

+ 5m 6D M 2 * + I 6G2 M 2 + B 6D M 2 + A5m 6D M 2 * [I] - 5m 4b L 2 M 2 * 

- B 4b L 2 M 2 + L 4s M 4b - L 4s 5m 2 M 2 * - L 4s B 2 M 2 - 5m 4a{l *)5m 2 M 2 * 

- 5m 4a ^)B 2 M 2 * - B 4a{ ii)B 2 M 2 - M 4b(2 *)5m 2 L 2 - M 4b B 2 L 2 - M 4b B 2 L 2 , 

+ I 4c M 2 I 2 + 5m 2 B 2 M 2 * L 2 i + 5m 2 L 2 5rh 2 * M 2 * + B 2 L 2 M 2 *5rh 2 r + B 2 L 2 M 2 B 2 i 
+ (B 2 ) 2 M 2 L 2/ 

M 29 = AM 29 + 2M %E B 2 + 2M QE[2 , ) 5m 2 + 2L 6M M 2 + 2I 6E2 M 2 

- 2L M{vy) B 2 M 2 - 2L 4s(y) B 2 M 2 - 2M 4a{2 « ) 5m 2 B 2 - M 4a (B 2 ) 2 - M 4a(2 * 2 *)(5m 2 ) 2 
+ 2L 2ll (B 2 fM 2 

M 30 = AM 30 + AM 6A I 2 + 2M 6B B 2 + 2M m2 *)5m 2 + 5m 6A M 2 , 

+ 2I 6B2 M 2 + B 6A M 2 + A5m 6A M 2 *[I] - 25m 4b{l « ) 5m 2 M 2 * - 25m 4b{v) B 2 M 2 * 

- 2B 4b{v) B 2 M 2 - 2M 4b(2 * ) 5m 2 B 2 - M 4b (B 2 ) 2 - M 46(2 * 2 *)(<5m 2 ) 2 + 2I 4s M 2 I 2 
+ 25m 2 B 2 5m 2 ,*M 2 * + {B 2 ) 2 M 2 »bm 2 „ + {B 2 ) 2 M 2 B 2 „ 

M 31 = AM 31 + 2L 6H3 M 2 

M 32 = AM 32 + L m2 M 2 + L 6G3 M 2 
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M 33 = AM 33 + 2L QG3 M 2 

M 34 = AM34 + L mi M 2 + L 6C3 M 2 + L Ax M Aa - 2L Ax L 2 M 2 

M 35 = AM35 + M m L 2 + L 6E3 M 2 + L 6G4 M 2 - 2L Ax L 2 M 2 

M 36 = AM 36 + M 6G L 2 + L 6B3 M 2 + L 6G5 M 2 - L Ac L 2 M 2 
+ L Al M Aa - 3L Al L 2 M 2 - M Aa (L 2 ) 2 + 2(L 2 ) 3 M 2 

M 37 = AM 37 + 2L QG2 M 2 

M 38 = AM 38 + AM m I 2 + 6m 6H M 2 * + B m M 2 + A5m m M 2 * [I] 
+ 2L Ax M Ab - 2L A Jm 2 M 2 * - 2L Ax B 2 M 2 

M 39 = AM 39 + Mq G L 2 + L 6G1 M 2 + L 6C2 M 2 + L Ac M Aa 

- 3L Ac L 2 M 2 - L Al L 2 M 2 - M Aa {L 2 f + 2{L 2 fM 2 

M 40 = AM 40 + AM 6G I 2 + Mq G L 2 + 5m 6G M 2 * + B 6G M 2 

+ A5m 6G M 2 * [I] - 5m Aa L 2 M 2 * - B Aa L 2 M 2 + L Ac M Ab - L Ac 5m 2 M 2 * 

- L Ac B 2 M 2 + L Al M Ab - L Al 5m 2 M 2 * - L Ai B 2 M 2 - 2M Ab (L 2 f 
+ 25m 2 (L 2 ) 2 M 2 * + 2B 2 (L 2 ) 2 M 2 

M 4 i = AM41 + 2M 6E L 2 + 2L 6C1 M 2 + 5rh 4a M 4a{2 * } + B Aa M Aa 

- 2B Aa M 2 L 2 , - AL As L 2 M 2 + A<5m 4a /4a(2*) - 2M 4a (2.)<5m 2 L2 - 2M Aa B 2 L 2 
+ AM Aa I Ax + 4B 2 L 2 M 2 L 2 , 

M 42 = AM 42 + AM 6C I 2 + 2M 6B L 2 + 5m 6G M 2 * + B 6G M 2 

+ A5m &c M 2 *[I] + 8m Aa M Abi2 *) - 5m Aa 5m 2 *M 2 * - 25m Ab L 2 M 2 * + B Aa M Ab 

- B Aa M 2 *5m 2 , - B Aa M 2 B 2 i - 2B Ab L 2 M 2 + A5m 4a / 4f) (2.) - 2M Ab ( 2 *- ) 5m 2 L 2 

- 2M Ab B 2 L 2 + AM Aa I Al + AM 4a (/ 2 ) 2 + 25m 2 L 2 5m 2 * M 2 * + 2B 2 L 2 M 2 *5m 2 , 
+ 2B 2 L 2 M 2 B 2 , 
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M 43 = AM 43 + M m B 2 + M m(3 * ) 5m 2 + 2L 6E2 M 2 + h H3 M 2 

- 2L 4x ( 2 t)B 2 M 2 

M 44 = AM44 + M 6G B 2 + M 6G{3 * ) 5m 2 + L 6E1 M 2 + L 6B2 M 2 

+ hc?,M 2 + L 4s M 4a — 2L 4s L 2 M 2 — L 4c ( 2 ')B 2 M 2 — L 4 i{ 2 ^B 2 M 2 

- M 4a B 2 L 2 i + 2B 2 L 2 M 2 L 2 , 

M 45 = AM45 + AM 6E I 2 + M 6C B 2 + M 6C(3 *)(5m 2 + 5m 6E M 2 * 

+ hczM 2 + B 6E M 2 + A5m 6E M 2 * [I] + 2L 4s M 4b - 2L 4s 5m 2 M 2 , 

- 2L 4s B 2 M 2 - 5m 4a ( 2 *)5m 2 M 2 * - 5m 4a ( 2 >)B 2 M 2 * - B 4a ^ 2 r ) B 2 M 2 - 2M 4b B 2 L 2 > 
+ I 4x M 2 I 2 + 25m 2 B 2 M 2 *L 2 , + 2(B 2 ) 2 M 2 L 2 , 

M 46 = AM 46 + M 6E B 2 + M 6E[3 * ) 8m 2 + 2L 6B1 M 2 + h E3 M 2 

+ 5rh 4b M 4a ( 2 *) + B 4b M 4a - 2B 4b M 2 L 2 , + A5m 4fe / 4a ( 2 *) - 2L 4s ^)B 2 M 2 

- M 4a{2 ,)5m 2 5m 2 * - M 4a{2 ,)B 2 5m 2l - M 4a B 2 B 2 < + AM 4b I 4x + I 4x M 2 I 2 
+ 2B 2 M 2 B 2 i L 2 i 

M 47 = AM47 + AM 6B h + M 6B B 2 + M 6B( 3*)5m2 + 5m 6B M 2 * 

+ I 6B3 M 2 + B &B M 2 + A5m 6B M 2 * [I] + 5rh 4b M 4b{2 *) - 5rh 4b 5rh 2 * M 2 * 

+ B 4b M 4b - B 4b M 2 *5rh 2 i - B 4b M 2 B 2 , - 5m 4b ( 2 *)5m 2 M 2 * - 5m 4b ^)B 2 M 2 * 

- B 4b{2 ,)B 2 M 2 + A5m 4b I 4b{2 *) - M 4b(2 ,)5m 2 5m 2 , - M 4b{2 *)B 2 8m 2 , - M 4b B 2 B 2 , 
+ AM 4b l 4l + AM 4b {l 2 f + 2l 4l M 2 l 2 + <5m 2 (5m 2 *) 2 M 2 * + B 2 5m 2 *M 2 *5m 2 , 

+ B 2 M 2 *5m 2 ,B 2 , + B 2 M 2 {B 2 ,f + M 2 (7 2 ) 3 

2. Divergence Structure of Quantities of Sixth- and Lower-Orders 

Substitution of M, into expresses the latter in terms of finite quantity AMj. However, 
still contains numerous (divergent) quantities of lower orders. To obtain a finite result 
the end we need the information on the divergence structure of these quantities, which 
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are listed in the following. The suffix 6a, (a = A,..., H), in L Qa , Bq , <5m 6cf refers to the 
diagrams similar to those of Fig. before H2 insertion is made. 

Vertex renormalization constants of sixth-order 

Lqai = hAi + AL 6A i + 2<5m 2 L 4s (i*) + 2B 2 L 4s ( V) - 5m 2 {5m 2 L 2 ^ 

+ B2L2'*) — i?2(5m2-^2'* + B2L2") + Lqai 
LqA2 = hA2 + AL 6A2 + L 2 L 4s + ^2-^4i(l*) + B 2 L^{ii) 

- L 2 (5m2-L 2 . + B 2 L 2 >) + hsL 2 + L 6A2 

L&A3 = hA3 + AL 6A3 + 2(5m 2 L 4s(1 *) + B 2 L 4s{1/) ) - 5m 2 (5m2L 2 .t* 

+ B2L2'*) — i?2(5m2-L2'* + B2L2") + L 6 A3 

Lqbi = hBi + AL 6B1 + 5m 2 L 4s (2») + /3 2 -^4s(2') + Srh 4b L 2 * 

+ -B45L2' - 5m25rri2*L2* - B 2 {5rh2>L2* + B 2 >L 2 >) + l/2J 6B + L 6B1 

L&B2 = I&B2 + ALq B2 + <5m 2 -L 4 i(2*) + -E>2-^4«(2') + L is L 2 

- B2L21L2 + /2-^4s — 12(5^2-^2* + B2L2') + -^6/32 

-^6/33 — -^6/33 + AL 6 £3 + I/ 2 I/4; + L^L 2 — (L 2 ) 2 L2 
+ -^2-^4/ — I2L2L2 + /4«-^2 + -^6B3 

-^6C*i = ^6Ci + ALgcx + 2L 2 L As + 5m 4a L 2 * + B Aa L 2 ' 

- 2L 2 (5m 2 L2» + £ 2 L 2 ' ) + 1/2J 6C + vL 6C i 

-^6C2 = -^6C2 + AL 6 C2 + / L 2 -£ 4 i + L ic L 2 — (I/ 2 ) 2 I/2 
+ -^2-^4c — I2L2L2 + L G C2 

L&C3 = -^603 + AL 6< 73 + L 4x L2 + hL^x + Lqc3 
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Lqdi — Iqdi + AL 6 di + 5m 2 L4 C (!.) + _B 2 -^4c(i') + L 2 L is 

— L 2 (5m 2 L 2 * + B 2 L 2 >) + L 6D1 

L%D2 — + ^Lq£> 2 + L 2 L ic + L 2 L^i — (L 2 ) 2 L 2 + Ii C L 2 + L§di 

Lqd3 = ^6D3 + AL 6 D3 + 8m 2 L Ac (i*) + B 2 L 4c ^) + L 2 L± S 

— L 2 (6m 2 L 2 * + B 2 L 2 i) + L6D3 

-^6D4 = -^6D4 + AL 6 £> 4 + 5^2-^4x(l*) + ^2-^4x(l') + -^6D4 

-^6D5 — -^6F>5 + AL 6 £> 5 + <5m 2 L 4c (3») + _B 2 L 4c ( 3 /) + L is L 2 

— B 2 L 2 iL 2 + J 2 AL 4s + L6D5 

-^6Bi = -^6Ei + ALesi + 8m 2 L^ c { 2 t) + B 2 L^ c ( 2 i) + La lS L 2 

— B 2 L 2 iL 2 + J 2 AL 4s + L 6 £! 

-^6F2 = ^6F2 + ALq E2 + 5m 2 L ix ( 2 *) + B 2 L 4:X (2') + -^6F2 
-^6F3 = -^6F3 + AL6S3 + -^4:r-^2 + L 2 L^ X + -^6^3 

-^6Fi = -^6Fi + AZ^fi + L 2 L,4 C + L^ C L 2 — (L 2 ) L 2 + J 2 AL 4c + I/6F1 

-^6F2 — -^6F2 + AL 6 ^ 2 + L 2 L ix + L^p 2 

-^6F3 = -^6F3 + ^^6F3 + 2Z/ 2 Z/ 4c — (L 2 ) 2 L 2 + I/6F3 

^6G1 — ^6G1 + AL 6G i + -^2-^4c + -^4c-^2 — (L 2 ) 2 L 2 + I 2 AL 4c + L 6G1 
-^6G2 = -^6G2 + AL 6 G2 + -^6G2 
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L%ga = ha + AL 6G 4 + L 2 L ix + L 6G 4 

-^6G5 = ^6G5 + AL 6G5 + -L 2 -L 4c + I/ 4 ;Z/ 2 — (L 2 ) 2 L 2 + I 2 AL A i + L 6G5 

-^6m — -^m + AL eH1 + L 4x L 2 + I 2 AL 4x + L^hi 

Lqhs = hm + AL 6 ^ 3 + Lqh3 
Wave-function renormalization constants of sixth-order 

5 

Bqa = — ^2 ^ 6Aa + ABqa + 2(<5m 2 -B 4 &(i») + B 2 B ib ^) 

a=l 

- 5m 2 (5m 2 B 2 ** + B 2 B 2 ,*) - B 2 (5m 2 B 2 i* + B 2 B 2 i<) + 2I 4s B 2 + £? 6j4 

5 

-B6B = — ^ -^6Bq + A_B 6B + 5m 2 B ib ( 2 *) + B 2 B ib ( 2 /) 

a=l 

+ 5m ib B 2 * + B ib B 2 / — 5m 2 5m 2 *B 2 * — B 2 (5m 2 rB 2 * + B 2 /B 2 i) 

+ I 2 B ib - I 2 (5m 2 B 2 * + B 2 B 2 i) + I 4l B 2 + 5 2 * [J]A5m 46 - J 6B + £ 6B 

5 

Bqc = — -^6Ca + ABqc + 2L 2 B ib + bm ia B 2 * + B ia B 2 ^ 

a=l 

- 2L 2 (5m 2 B 2 * + S 2J B 2 + / 2 5 4 „ - 2I 2 L 2 B 2 + 5 2 * [J]A5m 4a - J 6C + 5 6C 

5 

Bqd = — ^2 ^ Da + ^B§d + ^2-B 4a (i») + B 2 B ia (i') + I/ 2 -B 4b + L 4s -B 2 

a=l 

- L 2 (Sm 2 B 2 * + B 2 B 2 i) — B 2 L 2 /B 2 + J 4c -B2 + B 2 AL 4s + Bqd 
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5 

Bqe = — ^ hEa + ABq E + Sm 2 B 4a ^ 2 *) + B 2 B ia (2>) 
a=l 

+ 2L 4s B 2 — 2B 2 L 2 iB 2 + l 4x B 2 + 2B 2 AL 4s + Bqe 

5 

Bqf = — ^ hFa + ABq F + 2L 2 B 4a + 2L 4c B 2 — 3(L 2 ) 2 B 2 + 2B 2 AL 4c + B 6F 

a=l 

5 

-B6G = — ^2 ^ 6Ga + ^-^6G + L 2 B Aa + Z/ 4c -B 2 

0=1 

+ L 4l B 2 - 2(L 2 ) 2 B 2 + B 2 (AL 4l + AL 4c ) + B %G 

5 

Bqh = — ^ hn a + ABq H + 2L 4x B 2 + 2B 2 AL 4x + B eH 

a=l 

Mass renormalization constants of sixth-order 

SrriQA = ASttlqa + 2(5m 2 5rh 4b ^^ + B 2 5m 4b ^^) — 5m 2 (5m 2 5m 2 ** + B 2 5m 2 i*) 

— B 2 (5m 2 5m 2 i* + B 2 5fh 2 u) + Sitiqa 

SrriQB = A5m 6B + 5m 2 5m 4b ( 2 *) + B 2 5m 4b ( 2 >) + 8rh 4b 8m 2 * + B 4b 8fh 2 > 

— 5m 2 5rh 2 *5rh 2 * — B 2 (5m 2 '5m 2 * + B 2 i5fn 2 i) + I 2 5rh 4b — I 2 (5m 2 5rh 2 * 
+ B 2 5m 2 /) + <5m 6B 

^ m 6C = A<5m 6 c + 21/2<5m 4 b + <5m 4a <5rri2* + B 4a 8m 2 i — 2L 2 (5m 2 Srh 2 * 
+ B 2 bm 2 i) + I 2 5m 4a + 5m 6C 

5m 6Z) = A5m 6D + 5m 2 <5m 4a (i») + 5 2 5m 4 a(i') + L 2 5rh 4b - L 2 (5m 2 8fh 2 * 
+ B 2 5m 2 >) + 5m 6D 

<5m 6 £ = A5m 6E + 5m 2 5fn 4a ( 2 *) + B 2 bm 4a{2 >^ + 5m 6£ 

5m 6 ^ = A<5m 6i ? + 2L 2 5rh 4a + Sm^p 
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5m 6G = A5m GG + L 2 5m 4a + <5m 6 G 
5m 6H = A5m 6H + 5rh m 

Renormalization Scheme for Sixth-order Moment 

a 6A = M 6A - 28m 2 M 4b{1 * ) - 2B 2 M 4b + 5m 2 (5m 2 M 2 ,* + B 2 M 2 .) 
+ B 2 (5m 2 M 2 » + B 2 M 2 ) 

a eB = M 6B - 5m 2 M 46 (2*) - B 2 M 4b - 5m 4b M 2 * - B 4b M 2 
+ 5m 2 (5m 2 *M 2 * + B 2 *M 2 ) + B 2 (5m 2 M 2 * + B 2 M 2 ) 

a 6C = M 6C - 2L 2 M 4b - 5m 4a M 2 * - B 4a M 2 + 2L 2 (5m 2 M 2 * + B 2 M 2 ) 

a 6D = M 6D - L 4s M 2 - 5m 2 M 4a (i*) - B 2 M 4a - L 2 M 4b 
+ B 2 L 2 M 2 + L 2 (5m 2 M 2 * + B 2 M 2 ) + L 2 *5m 2 M 2 

a 6E = M 6E - 2L 4s M 2 - 5m 2 M 4a{2 ,) - B 2 M 4a + 2L 2 ,5m 2 M 2 + 2L 2 B 2 M 2 

a 6F = Mqp - 2L 4c M 2 - 2L 2 M 4a + 3L 2 L 2 M 2 

a 6G = Mqg - L 4c M 2 - L 4l M 2 - L 2 M 4a + 2L 2 L 2 M 2 

a m = M m - 2L 4x M 2 

M 6A = AMqa + 2<5to 2 M 46(1 * ) + 2B 2 M 4b - 5m 2 (5m 2 M 2 „ + B 2 M 2 *) 

- B 2 (5m 2 M 2 * + B 2 M 2 ) + 2I 4s M 2 

M 6B = AMqb + 5m 2 M 4b{2 * ) + B 2 M 4b + 5m 4b M 2 * + B 4b M 2 

- 5m 2 5m 2 *M 2 * - B 2 (5m 2 ,M 2 * + B 2 ,M 2 ) + l 2 AM 4b + M 2 * [I] A5m 4b 
+ {I 4l + (/ 2 ) 2 }M 2 
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Mq C = AM 6C + 2L 2 M 4b + 5m 4a M 2 » + B 4a M 2 - 2L 2 (5m 2 M 2 , 
+ B 2 M 2 ) + I 2 M 4a - 2I 2 L 2 M 2 + M 2 * [I] A5m 4a 

M 6D = AM 6D + L 4s M 2 + (5m 2 M 4a (i*) + B 2 M 4a + L 2 M 4b 
- B 2 L 2 ,M 2 - L 2 (5m 2 M 2 * + B 2 M 2 ) + J 4c M 2 

M 6E = AMqe + 2L 4s M 2 + 5m 2 M 4a(2 * ) + B 2 M 4a - 2L 2 ,B 2 M 2 + J^Mj 

M 6F = A.\ /,;/• + 2L 4c M 2 + 2L 2 M 4a - 3L 2 L 2 M 2 

M 6G = AMqg + L 4c M 2 + /. ;,.!/, + L 2 M 4a - 2L 2 L 2 M 2 

M m = AM m + 2L 4x M 2 

Various Fourth-Order Quantities 

Sm 4a = 5m 4a + A5m 4a 

5m 4b = 5rh 4b + A5m 4b + 5m 2 5m 2 * + B 2 5m 2 t 

L 4c ((vy) = -^4 C (i') 

La s {(v)') — L 4s (ii) 

L 4x = L 4x + I 4x + AL 4x 

L 4x — I 4x + AL 4x 

L 4c = L 4c + I 4c + AL 4c + L 2 L 2 

L 4c = I 4c + AL 4c + L 2 L 2 
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B 4a = B 4a — + AB 4a + 2L 2 B 2 — 2/ 4c 

B 4a = —I 4x + AB 4a + 2L 2 B 2 — 2I 4c 

L 4 i = L 4i + I 4l + (J 2 ) 2 + AL 4i + L 2 L 2 

L 4l = I 4l + (J 2 ) 2 + AL 4l + L 2 L 2 

L 4s = L 4s + 7 4s + AL 4s + 5m 2 L 2 * + £? 2 L 2 ' 

-B46 = + AB 4b + 5m 2 B 2 * + B 2 B 2 i + J 2 -B 2 — 2/ 4s — 7 4 ; 

M 4a(1 .) = AM 4a(1 ,) + L 2 M 2 * + / 4a(1 * } 

-^4o(2*) = AM 4a ( 2 .) + 7 4a ( 2 *) 

M 46(r) = AM 4b(1 *) + (5m 2 M 2 »* + B 2 M 2 «) + I 4b{1 * } 

M 4b(2 *) = AM 4K2 .) + 5m 2 *M 2 * + 7 46(2 *) + (7 2 M 2 * + M 2 .[/]<5m 2 .) - 2M 2 *[/]/ 2 

M 4a = AM 4a + 2 /.,.!/, 

M 4b = AM 4b + 5 2 M 2 + 5m 2 M 2 , + J 2 M 2 

AM 4a( 2.) = AM 4s - 2AM 4a(1 ») - 2AM 46(1 *) - AM 46(2 .) 

AB 4 = AB 4a + AB 4b 

AM 4 = AM 4a + AM 4b 

AL 4 = AL 4x + 2AL 4c + 2AL 4s + AL 4l 
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A5rri4 = A5m 4a + ASm^ 

2-B 4a (i») + -04a(2*) = — 4(3L 4c (i*) + L 4c ( 2 *)) — 4(L 4a; (i*) + L 4x ^*)) 
2-B 4 fe(i*) + -846(2*) = — 4(3L 4s (x*) + L 4s ( 2 *)) — 4(L 4i (!*) + L 4 ^ 2 *)) 

Various Relations Among Second-Order Quantities 

i?2' = B2 

= 

Srhy = Srri2 — 5m 2 ' 
Srri2» = S1712 — dmy 
B2** — — — 2^2Z^2**t ~t~ -^2*t*) 
L 2 '* = -^2* — I2* + AL 2 * 

1^2'* — — B2* ~~ ^^-^2* — 2^2* 

5m2'* = 5^2* — 5rri2'* 
5rri2* = 5rri2* + 61712* 
5rri2* = h + A<5m 2 * 
B 2 = —L 2 
B 2 = -h + Afi 2 
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